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1.  EXECUTIVE  SUMMARY 

The  objective  of  this  task  was  to  study  the  origin  of  the  phenomena 
that  contribute  to  the  overall  noise  in  photonic  phased  array  systems.  In 
particular,  the  task  included  the  study  and  the  characterization  of  the 
semiconductor  laser  low  frequency  noise  and  its  1/f  behavior. 

The  approach  adopted  for  this  study  included  three  parts.  In  the 
first,  a  comprehensive  study  of  the  effect  of  phase  noise  from  lasers  and 
other  sources  on  photonic  phased  array  systems  was  considered.  This  study 
considered  the  influence  of  additive,  as  well  as  multiplicative  noise  in 
photonic  systems,  and  resulted  in  the  important  conclusion  that  the  effect 
of  uncorrelated  multiplicative  phase  noise  in  a  phased  array  system  is 
diminished  as  the  number  of  array  elements  is  increased.  The  study 
further  revealed  that,  by  contrast,  the  additive  phase  noise  of  the  array  does 
not  diminish  with  the  number  of  elements.  Thus  for  an  array  with  a  large 
number  of  elements,  the  overall  signal-to-noise  ratio  will  be  independent  of 
the  size  of  the  array.  Results  of  this  segments  of  the  study  pointed  to  the 
significance  of  including  both  the  additive  and  multiplicative  noise  sources 
in  system  designs  of  photonic  phased  array  antennas. 

In  the  second  segment,  the  influence  of  modulation  on  the  noise  of 
semiconductor  lasers  was  experimentally  determined.  In  particular,  the 
dependence  of  the  low-frequency  noiseon  injection-current  modulation  in 
external  cavity  semiconductor  lasers  was  examined.  This  noise  can  limit  the 
performance  of  photonic  sensors  and  phased  array  antenna  systems,  and 
can  reduce  the  sensitivity  of  photonic  and  fiberoptic  sensors. 
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The  approach  followed  for  the  study  of  the  1/f  noise  in  the  spectrum 
of  semiconductor  lasers  was  the  consequence  of  the  need  to  examine  the 
noise  for  a  single  mode  laser,  based  on  the  findings  of  the  second  part  of 
the  study.  This  approach  included  a  model  of  the  laser  as  a  noise-driven 
resonant  amplifier.  A  dynamic  partial  wave  model  devised  for  this  system 
was  then  used  to  obtain  the  spectrum  of  the  output  field  of  the  laser.  This 
model  produced  interesting  and  suggestive  results,  including  a  power- 
independent  term  in  the  Schawlow-Townes  formula  for  the  linewidth  of  the 
semiconductor  laser,  resulting  from  the  multiplicative  noise  due  to  electron 
density  fluctuations  in  the  gain  medium.  The  model,  however,  failed  to 
directly  account  for  the  1/f  behavior  in  the  spectrum  as  a  consequence  of 
the  laser  action. 
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2.  INTRODUCTION 

Noise  is  typically  the  limiting  factor  in  the  performance  of  photonic 
RF  systems.  The  influence  of  noise  in  the  stability  of  fiber  optic 
distribution  links,  previously  studied  in  detail  at  JPL,  pointed  to  the  need 
for  a  comp-rehensive  examination  of  the  influence  of  the  phase  noise  of 
semiconductor  lasers  in  photonic  multi-element  arrays.  The 
Semiconductor  Laser  Low  Frequency  Noise  Characterization  task  was 
aimed  at  addressing  this  problem. 

The  approach  used  in  this  task  was  devised  to  systematically  consider 
the  influence  of  low  frequency  noise  of  semiconductor  lasers,  both 
experimentally  and  analytically,  and  then  develop  a  model  for  the  origin  of 
the  1/f  noise  in  the  laser  spectrum.  In  this  final  report  results  of  the  study 
will  be  presented.  The  presentation  is  divided  in  three  segments.  In  the 
first  part  the  effect  of  the  phase  noise  on  the  performance  of  phased  arrays 
is  analyzed.  In  the  second  part,  the  influence  of  modulation  on  the  laser 
noise  is  reviewed,  and  experimental  results  are  presented.  The  third 
segment  includes  the  presentation  of  a  model  for  the  noise  of  the 
semiconductor  lasers  which  yields  the  influence  of  the  multiplicative  noise 
on  the  linewidth  manifested  as  an  additional  term  to  the  Schawlow-Townes 
formula. 

While  the  study  did  not  reveal  the  origin  of  1/f  noise  in  the  spectrum 
of  semiconductor  lasers,  it  laid  important  groundwork  for  future 
examination  of  this  topic,  which  is  difficult  to  analytically  track,  yet  is  of 
fundamental  importance  in  the  performance  of  photonic  based  RF  systems. 
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Phase  Noise  in  Photonic  rf  Arrays 


The  function  of  rf  photonic  phased  arrays  is  fundamentally  based  on  the 
phase  control  of  individual  elements.  Thus  any  noise  due  to  individual  com¬ 
ponents  in  the  system  limit  the  system  performance.  Basically  the  total 
noise  in  the  system  is  determined  by  the  combination  of  both  additive  and 
multiplicative  components.  The  additive  term  is  produced  by  a  variety  of 
sources  including  the  thermal  noise  of  dissipative  elements,  the  shot  noise  of 
the  photodetectors,  and  the  laser  relative  intensity  noise  (RIN).  The  additive 
term  is  readily  measured  and,  in  most  instances,  is  easy  to  calculate;  so  its 
influence  on  the  performance  of  the  system  is  generally  straightforward.  The 
additive  noise  is  always  present  in  a  system,  and  can  be  determined  indepen¬ 
dent  of  the  presence  of  any  signal.  The  multiplicative  noise  originates  from 
low  frequency  gain  and  path  length  instabilities,  and  is  only  present  when 
the  signal  propagates  in  the  system.  This  term  usually  has  a  l/f“  behavior 
with  frequency,  and  thus  is  typically  dominant  in  frequency  regimes  close  to 
the-  carrier. 

Our  model  for  the  additive  and  multiplicative  noise  terms  deals  with  the 
uncorrelated  components.  The  model  is  given  in  detail  in  Appendix  I,  and 
essentially  generalizes  the  case  of  noise  for  a  single  element  with  unity  gain 
to  the  case  of  a  phased  array  with  M  elements.  The  analysis  leads  to  the 
conclusion  that  for  uncorrelated  white  Gaussian  noise,  the  SNR  for  an  M 
element  array  is  given  by. 
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where  SNR  is  the  signal-to-noise  ratio  of  the  M  element  array,  Eq  is  the 
amplitude  of  the  input  signal,  6(l){t)  is  the  multiplicative  phase  noise,  and 
N(t)  is  the  additive  noise  term.  In  the  above  expression  ()  represent  time 
averaging. 

From  the  above  equation  it  is  readily  concluded  that  the  additive  term 
in  the  SNR  remains  constant,  independent  of  the  number  of  elements  in  the 
array.  The  multiplicative  term,  by  contrast,  does  depend  on  the  number  of 
arrays,  and  is  reduced  by  factor  M.  As  M  increases,  the  multiplicative  noise 
term  reduces  to  a  negligible  factor,  and  the  noise  performance  of  the  array 
is  determined  by  the  additive  noise  contribution,  only. 

As  mentioned  above,  laser  RIN  is  a  multiplicative  noise  source.  Thus 
in  designing  photonic  rf  arrays  the  use  of  the  above  equation  permits  the 
determination  of  the  type  of  the  laser  best  suited  for  the  particular  system. 
Since  the  RIN  for  the  injection  current  modulated  lasers  is  generally  much 
larger  than  that  of  the  solid  state  YAG  lasers,  the  combination  of  the  additive 
noise  level  and  the  number  of  elements  required  will  dictate  which  type  will 


produce  the  highest  performance  in  a  given  array. 


4  Influence  of  Modulation  on  Semiconduc¬ 
tor  Laser  Noise 

In  this  task,  the  dependence  of  the  low  frequency  intensity  noise  on  injection- 
current  modulation  was  studied.  The  results  of  an  extensive  literature  search 
were  employed  to  determine  the  most  advantageous  manner  in  which  to  pro¬ 
ceed.  It  has  been  shown  previously  that  low  frequency  intensity  noise  (0  - 
1  GHz)  in  injection  modulated  semiconductor  lasers  is  upconverted  to  the 
vicinity  of  an  RF  modulation  signal.  The  heart  of  the  problem,  therefore,  lies 
in  determining  the  origin  of  the  low  frequency  fluctuations  in  the  laser  diode, 
and  how  the  low  frequency  noise  is  affected  by  injection  current  modulation. 

Two  competing  views  exist  regarding  the  origin  of  the  low  frequency 
fluctuations  in  semiconductor  lasers.  One  camp  contends  that  low  frequency 
intensity  fluctuations  are  due  to  trapping  of  carriers  in  the  semiconductor 
medium  with  a  1/f  power  law  frequency  spectrum,  while  the  other  view  holds 
that  competition  between  longitudinal  modes  of  the  laser  diode  causes  an 
enhancement  of  low  frequency  fluctuations.  The  initial  approach  we  adopted 
was  to  control  the  amount  of  coupling  between  longitudinal  modes  of  the 
laser,  and  to  observe  the  effect  on  the  low-frequency  fluctuations  and  close- 
to-carrier  RF  phase  noise. 

To  investigate  the  effects  of  mode  coupling,  an  external  cavity  semicon¬ 
ductor  laser  diode  was  constructed.  The  external  cavity  has  the  desirable 
effect  of  reducing  the  longitudinal  mode  frequency  separation  from  approxi¬ 
mately  450  GHz  to  1  GHz,  within  the  range  of  injection  current  modulation 
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frequencies.  Then,  the  mode-coupling  can  be  affected  by  injection  current 
modulation  at  appropriate  frequencies  corresponding  to  the  round  trip  light 
time  of  the  cavity  and  its  harmonics. 

Testing  of  the  external  cavity  semiconductor  laser  was  performed  under 
various  conditions  of  bias  current  and  modulation.  Preliminary  results  em¬ 
ploying  mode-locking  with  a  strong  sinusoidal  modulation  at  1  GHz  caused 
a  reduction  in  the  low  frequency  intensity  noise  of  the  laser  diode  in  ex¬ 
ternal  cavity  by  approximately  10  dB.  The  cause  of  the  noise  reduction  is 
presumed  to  be  the  reduction  of  competition  between  the  lasing  modes  of 
the  laser  when  mode-locked.  However,  an  unexpected  result  was  obtained 
when  the  1/2  sub-harmonic  of  the  mode-locking  frequency  (i.e.,  500  MHz) 
was  employed:  the  low  frequency  intensity  noise  was  reduced  by  another  10 
dB  for  frequencies  above  5  kHz.  These  results  were  presented  at  the  Third 
Annual  DARPA  Symposium  on  Photonics  Systems  for  Antenna  Applications 
(PSAA-III),  on  January  20th  1994  in  Monterey  California,  in  a  paper  enti¬ 
tled:  ’’Influence  of  Modulation  on  Noise  in  Semiconductor  Lasers,”  by  R.T. 
Logan  Jr.  and  L.  Maleki. 

5  Multiplicative  Noise  and  Laser  Linewidth 

A  theoretical  model  for  the  mode-competition  noise  reduction  in  a  multi- 
mode  laser  was  developed  and  computer  simulations  were  performed.  The 
model  is  based  on  a  Fabry-Perot  cavity  containing  a  medium  with  a  time- 
varying  index  of  refraction.  In  the  simulation,  the  center  frequency  and 
amplitude  of  the  dominant  mode  was  tracked  to  estimate  the  frequency  and 
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intensity  fluctuation  power  spectra.  The  computer  random-number  genera¬ 
tor  was  used  to  simulate  index  of  refraction  variations  with  a  white  power- 
spectral-density;  the  model  appeared  to  generate  l/f“  -type  noise  similar  to 
that  observed  in  the  laser,  and  proportional  to  the  noise  amplitude. 

This  initial  theoretical  result  served  as  tentative  confirmation  of  the  mode- 
competition  hypothesis.  Theoretically,  then,  the  1/f  noise  should  disappear 
completely  if  mode-  competition  was  completely  suppressed  by  perfect  mode¬ 
locking.  It  should  also  disappear  in  the  case  of  a  purely  single-mode  laser. 
However,  further  searching  of  the  literature,  and  experience  revealed  that  a 
residual  level  of  1/f  noise  is  often  observed  even  in  single-mode  distributed- 
feedback  laser  diodes.  The  computer  model  was  applied  to  the  single-mode 
laser  diode  again,  incorporating  index  of  refraction  variations  with  a  white 
power-spectral-  density.  Again,  1/f-type  frequency  and  intensity  fluctuations 
were  observed  in  the  simulated  output. 

To  verify  the  computer  result,  an  analytic  solution  was  derived  for  the 
problem  of  the  fluctuations  of  the  transmission  of  a  resonant  cavity  due  to 
variations  in  the  feedback  parameters.  The  analytic  model  is  an  extension 
of  a  linear  feedback  system  model  that  includes  the  effects  of  time- variation 
of  the  feedback  components.  In  this  linear  model,  the  transmitted  wave  at 
any  time  instant  is  comprised  of  a  summation  of  ’’partial  waves”  that  have 
existed  in  the  cavity  for  various  lengths  of  time.  In  the  case  of  no  feedback 
fluctuations,  the  partial  waves  may  be  summed  analytically,  and  the  familiar 
Fabry-Perot  cavity  result  is  obtained.  As  the  Q  of  the  cavity  is  increased, 
more  partial  waves  are  stored  in  the  cavity.  However,  when  the  feedback 
parameter  is  fluctuating,  the  output  wave  is  a  summation  of  many  partial 


waves  that  have  sampled  the  feedback  fluctuations  over  increasing  lengths  of 
time.  It  is  seen  through  this  analysis  that  in  a  high-Q  cavity,  the  fluctuations 
of  the  cavity  are  amplified  by  the  fact  that  many  partial  waves  make  up  the 
output  wave  at  any  instant,  and  each  partial  wave  carries  a  history  of  the 
cavity  fluctuations  with  it. 

The  analytic  model  results  serve  to  confirm  and  explain  the  computer  sim¬ 
ulation  results.  The  fluctuations  of  the  output  wave  due  to  refractive  index 
variations  with  white  power  spectral  density  acquire  a  1/fr  power  spectral 
density  at  frequencies  above  the  -3  dB  point  of  the  cavity,  but  remain  white 
at  frequencies  within  the  cavity  bandwidth.  This  noise  enhancement  effect 
may  explain  why  high-Q  cavities  typically  have  relatively  poor  long-  term 
stability. 

The  linear  model  is  not  sufficient  to  model  an  oscillator,  such  as  a  laser. 
However,  it  provides  interesting  new  insights  into  the  effects  of  noise  in  res¬ 
onators  and  resonant  amplifiers.  In  the  next  phase  of  work,  the  time- varying 
feedback  and  partial  wave  model  was  incorporated  into  a  non-linear  analysis 
of  the  resonant  cavity  above  the  oscillation  threshold.  This  formalism  en¬ 
tailed  modeling  the  laser  as  a  noise  driven  resonant  optical  amplifier  that  has 
random  gain  and  phase  fluctuations.  The  model  yields  an  additional  power 
independent  term  in  the  Schawlow-Townes  expression  for  the  linewidth  of 
the  laser.  As  the  power  is  increased,  the  model  predicts  a  re-broadening  of 
the  linewidth,  as  generally  observed  in  the  case  of  semiconductor  lasers. 

The  study  also  included  an  analysis  of  electron  density  fluctuations  in 
the  laser  gain  medium,  from  which  the  spectral  power  density  of  noise  due 
to  these  fluctuations  was  obtained.  Finally,  the  results  were  used  to  obtain 
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an  estimate  of  the  minimum  linewidth  of  the  semiconductor  laser. 

This  segment  of  the  work  is  the  basis  for  the  dissertation  of  R.  T.  Logan, 
which  will  be  submitted  in  the  future  in  partial  fulfillment  for  the  Ph.  D.  de¬ 
gree  in  Electrical  Engineering,  Electro-Physics  at  the  University  of  Southern 
California.  The  portion  of  the  dissertation  supported  in  the  task  is  presented 
in  Appendix  II. 
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Proceedings  of  SPIE  0/E  Ease  '94,  Optoelectronic  Signal  Processing  for  Phased- 
Array  Antennas  IV,  26-27  January  1994,  Los  Angeles,  CA. 

EFFECTS  OF  PHASE  NOISE 
FROM  LASERS  AND  OTHER  SOURCES 
ON  PHOTONIC  RF  PHASED-ARRAYS 

Ronald  T.  Logan  Jr.  and  Lute  Maleki 

Jet  Propulsion  Laboratory,  California  Institute  of  Technology 
4800  Oak  Grove  Drive,  Pasadena,  California  91109 

ABSTRACT 

The  beam  pattern  of  a  linear  phased-array  anterma  system  employing  a  photonic 
feed  network  is  analyzed  using  a  model  for  the  individual  feed  element  noise 
including  both  additive  and  multiplicative  equivalent  noise  generators.  It  is 
shown  that  uncorrelated  multiplicative  noise  power  of  the  individual  feeds  is 
reduced  by  a  factor  of  N  in  the  output  of  an  N-element  linear  array.  However,  the 
uncorrelated  additive  noise  of  the  individual  feed  paths  is  not  mitigated,  and 
therefore  will  determine  the  minimum  noise  floor  of  a  large  phased-array 
antenna. 

1.  INTRODUCTION 

In  phased-array  antennas,  the  beam  pattern  depends  critically  on  the  phase 
control  of  the  signals  at  the  individual  antenna  elements.  The  ability  to  feed  and 
adjust  the  phase  of  the  microwave  signals  to  the  individual  radiating  elements 
using  optical  fiber  and  photonic  components  offers -obvious  advantages  in  size, 
weight,  mechanical  flexibility,  and  cross-talk,  compared  to  metallic  waveguides 
and  phase-shifters.  Various  phased-array  antenna  system  architectures  with 
photonic  feed  networks  have  been  proposed,  however,  the  issues  of  phase 
stability  and  signal  purity  are  not  typically  addressed  in  these  proposals. 
Therefore,  determination  of  the  acceptable  phase  noise  contribution  of  the 
individual  active  feed  components  has  been  problematic.  In  this  paper,  the 
general  factors  contributing  to  the  phase  stability  of  an  array  feed  network  are 
outlined,  with  particular  attention  paid  to  the  type  of  noise  encountered  in 
photonic  feed  elements.  It  is  shown  that  the  analysis  of  array  phase  stability  must 
consider  both  additive  and  multiplicative  noise  generation  processes,  and  that 
the  additive  noise  of  the  active  feed  components  will  limit  the  phase  stability  of  a 
large  phased-array. 

2.  PHASED-ARRAY  SYSTEM  MODEL 

The  general  architecture  of  a  phased-array  antenna  system  comprised  of  M 
elements  is  depicted  schematically  in  Figure  1.  In  this  analysis,  the  phase  noise 
contribution  due  to  the  array  feed  and  antenna  elements  only  is  calculated.  The 
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effects  of  the  source  phase  noise  will  be  common  to  all  elements,  and  may  be 
treated  in  the  usual  manner  for  a  single  antenna  element. 


N„:  Additive  Noise  Generators 
&))„;  Multiplicative  Noise  Generators 


Figure  1.  Noise  Sources  in  M-Element  Phased- Array  Feed  System 


The  feed  system  is  driven  by  a  common  source  oscillator  whose  power  output  is 
divided  M  ways.  The  phase  delays  required  to  point  the  antenna  beam  are 
generated  in  the  separate  feeds.  To  study  the  effect  of  noise  on  the  elements,  we 
assume  a  simple  case  in  this  model:  the  signals  acquire  equal  delays  in  the  feeds, 
and  are  then  recombined  in  a  second  M-way  power  combiner.  Thus,  the  output 
signal  amplitude  is  scaled  to  be  equal  to  the  input  signal  amplitude,  to  facilitate 
comparison  between  a  single-element  antenna  system  and  an  array. 

3.  NOISE  PROPERTIES  OF  PHOTONIC  FEEDS  _ 


The  photonic  feed  elements  contain  active  components  such  as  laser  diodes, 
photodiodes,  and  amplifiers.  The  phase  noise  contribution  of  a  microwave  fiber 
optic  feed  system  is  therefore  comprised  of  an  additive  noise  term  and  a 
multiplicative  noise  term.  Laser  relative  intensity  noise  (RIN),  shot  noise,  and 
thermal  (Johnson)  noise  are  additive  noise  sources  which  are  present  at  all  times 
independent  of  signal  level.  Low  frequency  gain  or  path-length  instabilities  that 
modulate  the  microwave  signal  amplitude  and  phase  are  multiplicative  noise 
sources  that  are  only  observed  when  a  signal  is  present.  As  shown  in  Figure  2, 
additive  noise  usually  determines  the  noise  floor  at  higher  offsets  from  the  carrier 

frequency,  but  multiplicative  noise  often  has  a  1/f®  power  law  characteristic, 
l<a<2,  so  is  typically  dominant  close  to  the  microwave  carrier  frequency  [1]. 


A-3 


/ 


Figure  2.  Typical  Phase  Noise  of  Photonic  Feeci  System  at  10  GHz 


Additive  noise  sources  due  to  thermal  and  shot  processes  or  laser  RIN  are 
independent  random  processes  and  therefore  can  be  assumed  uncorrelated 
between  the  feed  elements.  Multiplicative  noise  may  or  may  not  be  uncorrelated 
between  elements,  depending  on  its  origin.  For  example,  thermal  expansion  or 
vibration  of  all  the  optical  fibers  in  the  feed  network  may  produce  a  phase 
modulation  that  is  common  to  all  elements,  whereas,  laser  or  amplifier-induced 
1/f  gain  fluctuations  will  be  uncorrelated  between  elements.  Noise  that  is 
common  to  all  the  elements  may  be  referred  to  the  source  oscillator  and  treated 
as  if  the  array  were  a  single  element.  In  the  analysis  that  follows,  all  of  the  noise 
sources  in  the  individual  feed  elements  are  assumed  to  arise  from  independent 
random  processes,  so  can  be  treated  as  uncorrelated  sources. 

It  is  noted  that  the  multiplicative  noise  is  not  detectable  by  a  standard  noise 
figure  measurement  in  the  microwave  signal  frequency  band.  In  fact,  it  is 
difficult  in  practice  to  predict  the  amount  of  multiplicative  noise  in  an  amplifier 
or  laser  diode,  because  the  noise  l^el  may  itself  be  a  function  of  the  modulation 
signal  frequency  or  amplitude.  Therefore,  the  amount  of  multiplicative  phase 
noise  is  usually  determined  empirically. 

4.  ANALYSIS 

We  now  proceed  to  calculate  the  effect  of  the  noise  added  by  the  feed  elements 
on  the  total  array  performance.  Consider  first  a  unity  gain  single  feed  element 
consisting  of  a  fiber  optic  link  and  electronic  amplifiers,  with  a  microwave  input 
signal  =  of  constant  amplitude  E„  at  microwave  frequency  co^.  The 
output  signal  amplitude  for  a  single  feed  element  may  be  written 

+  N'it)  (1) 

where  N'it)  represents  an  additive  Gaussian  noise  process  with  a  white  power 
spectral  density  from  dc  to  well  above  (W„,  and  5<pit)  is  a  small  multiplicative 
phase  noise  term.  The  multiplicative  noise  term  can,  in  general,  be  complex  and 
thus  also  represent  gain  fluctuations.  Factoring  out  the  sinusoidal  variation  at  co^ 
yields  the  slow  time-variation  of  the  output  field  around  the  microwave  carrier 
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=  E/^“’  +  N(f) 


(2) 


Limiting  the  analysis  to  a  band  of  frequencies  5co  in  width  around  co^,  the 
additive  noise  term  can  be  written  as  a  random  phasor: 

=  =  5rit)  + jSKt),  where  |5i(f)|  =  |5r(f)|,  and  are  assumed  to  be 

independent  Gaussian  noise  processes  with  white  power  spectral  density  from 
dc  to  dcoll.  Note  that  as  the  input  signal  amplitude  is  decreased  to  zero  in 
Equation  (2),  the  multiplicative  noise  term  vaiushes,  but  the  additive  noise  term 
N(t)  is  unchanged. 


The  total  array  noise  is  now  calculated  by  using  the  M  individual  feed  element 
expressions  from  Equation  (2)  in  the  standard  calculation  [2]  of  the  array  output 
field.  For  a  linear  array  at  a  steering  angle  (p,  the  field  distribution  in  the  far  field 
of  the  array  as  a  function  of  observation  angle  6  may  be  expressed  as 


E„,(0, «,  ()  =  +  N.(f) 


(3) 


Figure  3  is  the  calculated  radiation  pattern  for  a  linear  array  of  ten  antennas  vs. 
observation  angle  6  when  the  steering  angle  <p  =  0. 


Figure  3.  Calculated  radiation  pattern  vs.  observation  angle  for  a  10-element  linear  array. 

We  wish  to  investigate  the  magnitude  of  the  amplitude  and  phase  fluctuations  of 
the  main  lobe  versus  the  number  of  elements  in  the  array.  If  the  noise  sources  in 
the  feeds  are  uncorrelated,  the  statistics  of  the  time-variation  of  the  amplitude 
and  phase  of  the  main-lobe  peak  will  be  independent  of  the  steering  angle  (p. 
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Therefore,  for  the  purposes  of  this  noise  analysis,  the  array  can  be  modeled 
simply  as  M  equal-length  feed  systems  sandwiched  between  two  back-to-back 
ideal  M-way  power  splitters,  as  depicted  in  Figure  1.  Now,  the  output  field 
amplitude  is  equal  to  the  input  field  amplitude,  independent  of  the  number  of 
elements  M. 


At  the  main  lobe  peak,  ;rcos  6  =  <f>,  and  the  signal  amplitude  at  the  output  of  the 
M-way  coupler  in  Figure  1  corresponds  to  the  main  lobe  peak  amplitude  of  the 

antenna  pattern,  divided  by  Vm.  For  small-angle  phase  noise  S(p{t) « 1  radian, 
the  time  variation  of  the  output  field  can  be  written 


First  consider  the  output  of  an  "array"  comprised  of  only  one  feed  element  with 
unity  gain.  In  the  above  equation,  this  corresponds  to  the  case  of  M=l.  The 
output  field  is  then  given  by 

E^,it)  =  E,  +  j5<l>(t)E,+N{t).  (5) 

which  is  just  the  same  as  the  expression  for  a  single  feed  given  earlier  in  Equation 
(2),  as  required.  The  signal  and  noise  power  are  proportional  to  the  time-average 
of  the  squared-magnitude  of  the  output  voltage: 

r(E,£.-)  +  (£,(;WE;)+(EX) 

“  +{is^Ex)+{m,(iSfyE:)+{is^E,N-) 

+{ne;) + (n(;5«)'e;) + (nn-) 

-where  the  angle-brackets  (  )  denote  the  time-avorage  of  the  enclosed  quantity, 
and  the  explicit  time-dependence  of  the  random  functions  has  been  dropped  for 
clarity.  For  independent  zero-mean  noise  processes,  the  time-average  of  products 
of  the  constant  and  random  terms  are  zero,  because  we  have  assumed  zero-mean 
random  noise  processes.  Also,  the  average  of  the  product  of  two  uncorrelated 

terms,  such  as  is  zero.  But  the  time-average  of  the  square  of  any  signal 

or  noise  term  is  non-zero.  Thus,  the  ratio  of  signal  to  noise  power  (SNR)  for  one 
feed  is 


c\jf?\  — _ V  / _  (y) 

Similarly,  for  the  case  of  M  parallel  feeds  with  independent  equal-amplitude 
multiplicative  and  additive  noise  sources,  the  output  field  is  given  by 
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The  individual  powers  of  equal-amplitude,  independent  (therefore  uncorrelated) 
noise  sources  may  be  added  linearly.  Thus,  the  equivalent  noise  voltage  due  to 

the  sum  of  M  uncorrelated  equal-amplitude  noise  sources  is  just  ^^M  times  the 
amplitude  of  a  single  noise  source.  Since  all  cross-terms  between  uncorrelated 
noise  sources  average  to  zero,  we  can  write  the  output  field  for  an  M-element 
array  in  terms  of  a  single  multiplicative  noise  source  d(p(.t)  and  a  single  additive 
noise  source  N(t) 


£«(t)  =  £.  +  ^’^^  +  N(f).  (9) 

Now,  the  ratio  of  signal  to  noise  power  in  the  combined  output  of  an  M-element 
array  is 


(lW»f|£.f)  ^ 

M 


(10) 


The  multiplicative  noise  power  is  thus  mitigated  by  a  factor  of  M  in  the 
combined  output  of  an  M-element  array,  whereas  the  signal  power  and  additive 
noise  power  are  unchanged  from  the  single-element  case.  Figure  4  illustrates 
these  results  for  an  array  of  ten  elements  compared  to  a  single  element. 


LOG(OFFSET  FREQUENCY) 

Figure  4.  Comparison  of  Single-element  and  10-element  Phase  Noise 


The  close-to-carrier  1/f®  multiplicative  noise  of  the  combined  array  output  is 
reduced  by  10  dB  compared  to  the  noise  of  an  individual  element.  The  white 
additive  noise  power  further  from  the  carrier  is  unchanged  from  the  single¬ 
element  case.  This  behavior  is  analogous  to  the  improved  frequency  stability 
obtained  from  an  ensemble  of  oscillators,  compared  to  the  stability  of  a  single 
oscillator.  This  property  may  therefore  make  phased-array  antennas  more 
desirable  than  single-element  anteimas  for  applications  in  which  high  levels  of 
long-term  phase  stability  are  required.  Alternatively,  this  property  relaxes  the 
requirements  on  multiplicative  phase  noise  for  the  elements  of  a  large  phase 
array.  It  is  emphasized  that  the  additive  noise  requirements  are  not  relaxed, 
however. 
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5.  SUMMARY 


It  was  shown  that  as  the  number  of  array  elements  M  is  increased,  the  effect  of 
uncorrelated  multiplicative  phase  noise  of  the  feed  elements  on  the  total  array 
stability  is  diminished.  However,  the  uncorrelated  additive  noise  of  the  feed 
elements  is  not  diminished,  so  that  the  signal-to-noise  ratio  becomes  independent 
of  array  size  for  large  enough  M.  It  is  therefore  important  to  quantify  both  the 
additive  and  multiplicative  noise  of  the  feed  elements  to  correctly  predict  the 
total  array  phase  stability. 
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Abstract 


A  general  formalism  is  developed  for  analyzing  the  output  field 
fiuctuations  of  a  resonant  optical  amplifier  that  has  random  gain  and 
phase  fluctuations  by  a  modification  of  the  partial  wave  model.  The 
partial  wave  model  provides  an  intuitive  physical  picture  for  the  effect 
of  multiplicative  noise  in  the  optical  resonant  amplifier.  Then,  by  mod¬ 
elling  a  laser  as  a  noise-driven  resonant  optical  amplifier,  it  is  shown 
that  multiplicative  noise  generates  an  additional  power-independent 
term  in  the  Schawlow-Townes  formula  for  the  linewidth  of  the  laser, 
and  a  re-broadening  of  the  linewidth  at  high  output  power  levels,  as  is 
typically  observed  in  single-mode  semiconductor  lasers  at  high  power 
levels.  Time-varying  complex  gain  constitutes  a  multiplicative  gain 
and  phase  noise  source  that  is  transformed  to  the  output  field  in  a 
different  way  fi:om  additive  noise  sources,  such  as  spontaneous  emis¬ 
sion.  The  formalism  developed  allows  calculation  of  the  transmitted 
field  power  spectra  of  amplitude  and  phase  fluctuations  due  to  multi¬ 
plicative  noise.  It  is  shown  that  the  effect  of  multiplicative  noise  on 
the  transmitted  optical  field  is  enhanced  as  the  net  round-trip  geiin  in 
the  resonator  is  increased,  so  that  the  multiplicative  noise  ultimately 
determines  the  minimum  linewidth  for  the  resonant  amplifier  in  the 
high-gain  limit.  A  det^uled  analysis  of  electron  density  fluctuations  in 
a  semiconductor  laser  gain  medium  is  also  performed,  from  which  the 
power  spectral  density  of  multiplicative  phase  noise  due  to  electron 
number-density  fluctuations  in  the  gain  medium  is  obtmned.  This  re¬ 
sult  is  then  used  in  the  multiplicative  laser  noise  model  to  obtain  an 
estimate  of  the  minimum  linewidth  of  a  semiconductor  laser.  Although 
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the  multiplicative  noise  analysis  is  applied  to  a  resonant  optical  ampli¬ 
fier,  the  formalism  is  general  and  should  be  applicable  to  the  analysis 
of  other  types  of  feedback  systems  perturbed  by  multiplicative  noise. 
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1  Introduction  and  Overview 


A  formalism  is  developed  for  analyzing  the  output  field  fluctuations  of  a 
resonant  optical  amplifier  with  a  time-varying  complex  gain  medium,  by  a 
modifled  partial-wave  model.  A  linear  model  of  a  laser  as  a  noise-driven  res¬ 
onant  optical  amplifier  is  then  modified  using  the  results  of  the  time-varying 
partial  wave  analysis.  This  new  laser  model  provides  an  intuitive  picture  of  a 
mechanism  for  power-independent  linewidth  and  linewidth  rebroadening  in 
semiconductor  lasers. 

In  a  resonant  amplifier,  time- varying  phase  and  gain  constitute  multi¬ 
plicative  noise  sources  whose  effects  are  transformed  to  the  output  field  in  a 
different  way  from  additive  noise  sources,  which  are  typically  treated  as  ad¬ 
ditional  inputs  to  the  system.  The  inclusion  of  multiplicative  gain  and  phase 
fluctuations  in  a  partial-wave  analysis  of  an  optical  resonant  amplifier  is  the 
principal  novelty  of  this  work.  The  partial  wave  model  provides  a  simple 
physical  picture  for  the  effect  of  multiplicative  noise  in  the  optical  resonant 
amplifier.  The  resonant  amplifier  is  modeled  as  a  delay-type  feedback  system 
(Figure  1)  with  input,  output,  and  net  round-trip  loop  gain  less  than  unity, 
as  shown  in  Figure  1.  The  output  field  at  any  time  instant  is  comprised  as 
a  sum  of  partial  waves;  the  “oldest”  partial  waves  have  travelled  more  times 
around  the  loop,  so  that  the  “memory”  of  the  cavity  extends  fiurther  back  in 
time.  When  no  multiplicative  noise  is  present,  the  width  of  the  cavity  reso¬ 
nance  peak  decreases  as  the  net  gain,  and  hence  the  cavity  memory  time,  is 
increased.  However,  in  the  presence  of  multiplicative  noise,  the  older  partial 
waves  accumulate  progressively  larger  random  phase  fluctuations,  leading  to 
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increasing  output  field  fiuctuations.  For  high  enough  net  round-trip  gain, 
the  multiplicative  noise-induced  fiuctuations  will  dominate  the  output  field 
fluctuations,  and  the  linewidth  will  not  derease  with  further  increases  in  net 
gain.  Eventually,  as  the  gain  is  increased  further,  the  later  partial  waves 
accumulate  so  much  phase  that  they  begin  to  destructively  interfere  with 
the  earlier  partial  waves,  leading  to  large  output  fluctuations,  and  increasing 
linewidth. 

By  viewing  the  laser  as  a  noise-driven  resonant  amplifier,  and  including 
the  results  of  the  partial- wave  analysis,  a  model  is  obtained  for  the  behavior 
of  the  laser  with  multiplicative  noise  as  the  net  round-trip  gain  is  varied.  It 
is  shown  that  multiplicative  noise  leads  to  an  additional  power-independent 
regime  in  the  Schawlow-Townes  formula  for  the  linewidth  of  a  laser  at  in¬ 
termediate  power  levels,  followed  by  linewidth  rebroadening. as  the  power  is 
increased  further.  From  this  result,  an  estimated  value  of  120  Hz  for  the 
miniTmim  power-independent  linewidth  of  a  typical  single-mode  DFB  semi¬ 
conductor  laser  due  to  electron  number-density  fluctuations  is  obtained. 

Multiplicative  noise  arises  in  an  optical  feedback  system  when  a  multi¬ 
plicative  factor  affecting  the  field,  such  as  gain  or  round-trip  phase,  fluctuates 
in  time.  Additive  noise  arises  from  photons  which  are  added  randomly  to 
the  field,  such  as  by  spontaneous  emission.  In  addition  to  the  fundamen¬ 
tal  addtive  noise  due  to  spontaneous  emission,  sources  of  multiplicative  gain 
and  phase  fluctuation  will  always  be  present  in  the  components  of  a  laser 
system  at  some  level.  It  has  been  shown  [12]  that  statistical  fluctuations  in 
temperature  and  density  can  have  significant  effect  on  the  optical  phase  of 
a  wave  passing  through  an  optical  fiber  at  a  constant  temperature  T.  For 


example,  in  a  fiber  laser,  this  source  of  fluctuation  constitutes  a  fundamental 
source  of  multiplicative  phase  noise  that  arises  from  the  same  thermodyamic 
considerations  responsible  for  the  unavoidable  Johnson  noise  in  electronic 
components,  and  is  a  basic  consequence  of  the  Fluctuation-Dissipation  The¬ 
orem  [16].  In  other  laser  systems,  the  net  gain  and  optical  length  of  the  laser 
cavity  can  fluctuate  due  to  dye-stream  width  fluctuations,  mirror  vibrations, 
pump  power  fluctuations,  or  other  technical  factors.  In  typical  laser  systems, 
it  is  these  technical  sources  of  multiplicative  noise  that  determine  the  actual 
observed  linewidth  [1],  and  not  fundamental  additive  quantum  noise  due  to 
spontaneous  emission,  as  predicted  by  the  well-known  modified  Schawlow- 
Townes  formula  [19],  [21]. 

The  Schawlow-Townes  formula  for  laser  linewidth  predicts  an  inverse  de¬ 
pendence  on  the  output  power  given  as 


6v  =  {dvff-K 


hu  No 


(1) 


P  N2-N1 

where  6uc  is  the  “cold  cavity”  linewidth  determined  by  the  losses,  h  is 
Planck’s  constant,  u  is  the  oscillation  frequency,  P  is  the  output  power, 
and  A^2,  Ni  are  the  populations  of  the  upper  and  lower  states  of  the  atomic 
transition  responsible  for  the  optical  gain.  For  large-scale  gas  lasers  with 
high-reflectivity  mirrors,  equation  (1)  predicts  linewidths  on  the  order  of 
several  hertz  for  output  powers  in  the  milliwatt  range.  In  real  lasers,  the 
observed  linewidth  typically  several  orders  of  magnitude  larger  than  this 
quantum  linewidth  prediction,  due  to  technical  sources  of  noise. 

A  notable  exception  is  the  semiconductor  laser.  The  combined  effects 
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of  large  output  coupling  (i.e.,  low  Q  cavity)  and  amplitude-phase  coupling 
through  the  carrier  density  [2],  [3]  yield  a  predicted  fundamental  quantum 
linewidth  in  the  kHz  to  100  MHz  range  for  typical  power  levels.  The  pre¬ 
diction  of  the  appropriately  modified  linewidth  equation  [2]  for  typical  facet 
reflectivity  of  32  percent  yields  an  estimated  linewidth  of  approximately  30 
MHz  at  a  power  level  of  10  milliwatts,  and  experiments  are  typically  in  good 
agreement  with  theory  for  lower  power  levels.  However,  a  power-independent 
linewidth  of  several  MHz  [11],  [22]  and  linewidth  rebroadening  [6],  [5]  in 
single-mode  semiconductor  lasers  is  typically  observed  at  high  output  power, 
which  is  not  predicted  by  the  modified  Schawlow-Townes  relation.  The  excess 
linewidth  is  detrimental  to  applications  requiring  a  high  degree  of  coherence, 
such  as  atomic  spectroscopy  and  coherent  communications  systems. 

The  origin  of  the  excess  linewidth  in  semiconductor  lasers  is  not  so  clear 
as  in  the  case  of  large-scale  lasers,  since  the  typical  macroscopic  technical 
sources  of  noise,  such  as  mirror  vibrations,  dye  stream  fluctuations,  and  the 
like  are  absent.  Early  work  on  the  power-independent  linewidth  attributed 
it  to  fundamental  multiplicative  refractive  index  fluctuations  due  to  electron 
number-density  fluctuations  [11].  In  that  work,  the  transformation  of  refrac¬ 
tive  index  variations  6n  to  fluctuations  of  the  laser  frequency  6i>  was  derived 
from  the  phenomenological  relation 


V  n 

The  refractive  index  fluctuations  were  included  in  a  root-mean-square  sense, 
derived  from  thermodynamic  considerations  regarding  the  number-density 
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fluctuations  of  the  electrons  contained  in  the  active  medium.  However,  this 
model  can  not  explain  the  onset  of  linewidth  rebroadening.  In  later  work, 
various  other  mechanisms  have  been  advanced  to  explain  the  origin  of  the 
power-independent  linewidth,  such  as  1//  noise  [22],  and  spatially-dependent 
temperature  and  carrier  fluctuations  [24].  Linewidth  rebroadening  has  been 
explained  in  terms  of  excess  carrier  density  in  the  confinement  regions  of 
quantum  well  lasers  [5],  and  spatial  hole  burning  [6].  These  or  other  mech¬ 
anisms,  such  as  injection-current  fluctuations,  may  be  operative  in  various 
combinations  in  any  particular  laser.  At  present,  however,  there  appears  to 
be  no  uniform  explanation  of  a  mechanism  responsible  for  both  the  power- 
independent  linewidth  and  linewidth  rebroadening  in  single-mode  semicon¬ 
ductor  lasers. 

The  dissertation  is  organized  as  follows:  In  Chapter  2,  a  discrete-time  for¬ 
mula  for  the  slowly-varying  complex  envelope  of  the  output  field  is  obtained, 
equation  (16).  Two  regimes  of  the  output  fluctuations  are  examined:  small- 
angle  and  large-angle.  The  small-angle  approximation  yields  equation  (30), 
which  has  the  form  of  a  finite-impulse-response  filter  acting  on  the  multi¬ 
plicative  and  input  noise  processes  to  produce  the  output  fluctuations.  Prom 
this  expression,  the  impulse  responses  and  frequency  responses  of  the  reso¬ 
nant  amplifier  for  multiplicative  and  input  fluctuations  are  obtained,  finally 
yielding  the  output  noise  power  spectrum,  equation  (46).  In  the  limit  of  high 
net  round-trip  gain,  this  becomes  equation  (47),  and  the  output  fluctuation 
power  spectrum  is  dominated  by  multiplicative  noise.  This  is  the  mechanism 
for  the  power-independent  linewidth.  Numerical  simulation  results  for  the 
large-angle  regime  illustrate  increasing  amplitude  and  phase  fluctuations  due 
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to  the  random- walking  phase  and  gain  of  the  partial  waves  Figures  2  - 
19.  With  increasing  net  gain,  the  probability  distribution  of  the  output  in¬ 
tensity  in  the  large  angle  regime  moves  to  lower  values  as  in  17  -  2  ^  and  the 
power  spectrum  broadens, (22  -  26  g'his  is  the  mechanism  for  the  linewidth 
rebroadening. 

In  Chapter  3,  for  the  small-angle  regime,  the  output  field  power  spectrum 
of  the  resonant  amplifier  is  found  to  have  a  Lorentzian  Uneshape,  that  has 
width  given  by  equation  (52)  for  a  multipUcative  noise  source  that  has  white 
power  spectrum.  This  result  has  the  same  form  as  the  phenomenological 
result  of  equation  (2).  In  the  large-angle  regime,  the  linewidth  is  predicted 
to  broaden,  which  is  not  predicted  by  equation  (2).  The  resonant  ampli¬ 
fier  analysis  is  finally  appUed  to  the  laser  by  modelling  it  as  a  noise-driven 
resonant  amplifier.  The  linewidth  of  the  laser  has  a  power-independent  min¬ 
imum  for  the  small-angle  regime,  and  rebroadens  following  the  onset  of  the 
large-angle  regime.  A  minimum  linewidth  for  a  typical  semiconductor  laser 
is  estimated  to  be  120  Hz  due  to  electron  density  fluctuations. 
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2  General  Analysis  of  Resonant  Feedback 
System  with  Multiplicative  Noise 

In  this  Chapter,  the  standard  linear  analysis  of  a  delay-type  feedback  sys¬ 
tem  is  generalized  to  treat  time- varying  gain  and  phase  perturbations  of  the 
components  in  the  feedback  system.  The  gain  and  phase  coefficients  multi¬ 
ply  the  complex  amplitude  of  the  field  as  it  passes  through  the  amplifying 
medium,  so  that  gain  and  phase  fluctuations  constitute  sources  of  multiplica¬ 
tive  noise.  First,  to  frame  the  discussion  and  define  terminology,  the  basic 
analysis  of  a  linear  delay-type  feedback  system  without  multiplicative  noise 
will  be  reviewed,  and  applied  to  the  case  of  a  resonant  optical  amplifier.  The 
gain  medium  is  treated  as  a  linear  amplifier  operating  below  saturation.  This 
amounts  to  a  review  of  the  standard  treatment  [4]  of  a  resonant  optical  ampli¬ 
fier,  from  which  the  familiar  Airy  formula  is  obtained.  Next,  a  discrete-time 
formalism,  valid  for  a  single  oscillation  mode  at  a  cavity  resonance  frequency, 
is  developed  including  multiplicative  phase  and  gain  fluctuations  of  the  am¬ 
plifying  medium  between  the  cavity  mirrors. 

In  the  standard  analysis  of  a  resonant  cavity,  the  transmitted  output 
field  is  computed  from  an  infinite  sum  of  partially-reflected  fields  inside  the 
cavity,  or  “partial  waves.”  Cavity  resonances  occur  at  discrete  frequencies 
when  all  of  the  partial  waves  add  constructively.  For  a  perfectly  stable  cavity 
with  net  round-trip  gain  less  than  unity,  an  analytic  form  for  this  summation 
obtains.  However,  in  the  presence  of  multiplicative  noise,  the  partial  waves  in 
the  resonator  will  have  randomly-varying  phase  and  amplitude,  so  the  usual 
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analytic  summation  is  not  applicable.  Therefore,  in  the  present  analysis,  the 
output  field  is  computed  as  a  random  phasor  summation  of  all  the  partial 
waves  contained  in  the  cavity  at  discrete  instants  of  time.  The  case  of  an 
ideal  noiseless  input  field  is  examined  first  to  appreciate  the  effect  of  the 
multiplicative  noise  on  the  transmitted  field.  Then  the  more  realistic  case  of 
an  input  field  with  phase  and  amplitude  noise  is  treated. 

Two  regimes  are  examined:  a  small-angle  linear  regime,  and  a  large- 
angle  regime.  In  the  small-angle  regime,  the  fiuctuating  complex  phase  of 
the  partial  waves  is  less  than  0.1  radians  about  the  mean.  In  this  regime, 
expressions  for  the  output  field  amplitude  and  phase  are  analytically  derived 
in  terms  of  the  power  spectra  of  the  multiplicative  gain  fluctuations  and  ad¬ 
ditive  input  fluctuations,  and  the  net  round-trip  gain.  But  it  is  seen  that 
the  multiplicative  noise  is  enhanced  by  the  round-trip  net  gain,  and  that  the 
resonant  amplifier  system  acts  like  a  linear  finite-impulse-response  filter  in 
transforming  the  multiplicative  fluctuations  and  the  input  field  fluctuations 
to  the  output  field.  The  impulse  responses  of  the  resonant  amplifier  for  both 
multiplicative  noise  and  input  field  noise  are  found  to  decay  exponentially, 
leading  to  Lorentzian  terms  in  the  output  field  power  spectrum.  The  fa¬ 
miliar  filtering  characteristics  of  the  Fabry-Perot  cavity  are  obtained  for  the 
input  additive  fluctuations.  The  response  for  the  multiplicative  fluctuations 
is  enhanced  as  the  gain  is  increased,  whereas  the  response  for  the  input  field 
fluctuations  decreases  with  increasing  gain.  For  large  enough  gain,  the  mul¬ 
tiplicative  noise  dominates  the  output  field  fluctuations.  It  is  shown  that 
the  enhancement  of  the  multiplicative  fluctuations  limits  the  output  power 
spectrum  to  a  minimum  width  determined  by  the  net  round-trip  gain  and  the 
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multiplicative  noise  strength,  in  contrast  to  the  prediction  of  the  standard 
analysis. 

For  the  case  of  large-amplitude  multiplicative  noise  exceeding  the  limits 
of  the  small-angle  approximation,  numerical  simulations  are  performed  using 
simulated  noise  time  series.  In  this  regime,  the  power  transmission  of  the 
cavity  with  multiplicative  phase  fluctuations  grows  progressively  more  erratic 
with  increasing  net  gain.  For  large  enough  net-gain,  the  phase  of  the  output 
field  “wraps  around”  more  than  27r  radians,  causing  discontinuous  jumps 
of  the  output  phase,  and  fast,  spiky  output  amplitude  fluctuations  on  time 
scales  much  faster  than  the  cavity  photon  lifetime.  Also,  the  probability 
distribution  of  the  output  intensity  becomes  more  like  that  of  a  thermal 
source,  and  the  power  spectrum  becomes  broadened.  In  this  regime,  the 
impulse  response  is  no  longer  exponentially  decaying,  but  becomes  erratic 
and  enhanced  at  large  times,  resulting  in  enhancement  of  the  low-frequency 
portion  of  the  output  field  power  spectrum. 

Finally,  it  is  noted  that  although  the  primary  concern  of  this  analysis  is 
the  optical  resonant  amplifier,  the  resuls  of  this  analysis  should  be  applicable 
to  noise  in  other  types  of  feedback  systems. 

2.1  Review  of  Standard  Linear  Feedback  Theory 

To  begin  the  development,  the  standard  analysis  of  a  simple  delay- type  feed¬ 
back  system  with  linear,  time-invariant  components  will  be  reviewed,  and 
then  applied  to  a  resonant  optical  amplifier  similar  to  the  approach  found  in 
Siegman  [4].  Consider  the  feedback  system  model  of  an  optical  cavity  illus- 
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trated  in  Figure  1.  The  input,  output,  and  internal  states  are  electric  fields 
with  sinusoidal  time  dependence.  The  input  field  R  incident  from  the  left  of 
Figure  1  is  assumed  to  be  a  sinusoid  of  frequency  u>  and  amplitude  Ar, 

R{lj)  =  Ah  (3) 

The  output  field  C{u)  exiting  to  the  right  also  has  the  same  sinusoidal  depen¬ 
dence.  It  is  desired  to  compute  the  transfer  characteristics  of  the  feedback 
system  from  the  input  state  R{u!)  to  the  output  state  C{uj). 

The  system  has  an  input  coupler  1  on  the  left  of  Figure  1,  and  an  output 
coupler  2  on  the  right.  These  couplers  have  branching  ratio  pi,2  for  the 
reflected  amplitude,  and  71,2  for  the  transmitted  amplitude.  The  couplers 
are  assumed  to  be  lossless,  so  that  712  +  Pi, 2  “  system  has  complex 

forward  gain  coefficient  G{u))  and  reverse  gain  coefficient  Hiui).  A  crucial 
assumption  is  that  the  gains  G[u})  and  H{uj)  are  linear,  so  that  the  output 
field  can  be  computed  as  a  vector  addition  of  partial  waves.  If  the  gain  is 
non-linear,  then  mixing  products  will  be  generated  at  new  firequencies.  We 
will  assume  that  the  cavity  behaves  as  a  linear  resonant  amplifier,  with  the 
output  field  amplitude  linearly  related  to  the  input  field  amplitude. 

Throughout  the  remainder  of  the  analysis,  the  explicit  sinusoidal  time 
dependence  of  R{uj),  G{uj)  etc.,  will  be  understood,  and  R,  G  etc.  will  be 
used  instead.  For  a  pure  sinusoidal  input  field  R  that  has  been  applied  for  an 
infinite  time,  the  steady-state  output  field  C  is  the  summation  of  an  infinite 
number  of  partial  waves  which  may  be  written 

G  =  (7i72)Gi? -b  (7i72)(piP2)G^-f^-R  +  (7i72)(PiP2)^G'^-f^^-R 
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Hii'y2){piP2YG^H^R  +  --- 

=  (7i72)Gi?[l  +  P1P2GH  +  {P1P2GHY  +  {P1P2GHY  +  •  •  •]  (4) 

The  term  outside  the  square  brackets  represents  the  transmission  of  the  input 
and  output  couplers  7171,  and  the  forward  gain  G  for  first  one-way  pass 
through  the  cavity.  The  terms  inside  the  brackets  represent  the  successive 
feedback  terms  which  experience  a  net  round-trip  gain  P1P2GH  on  each  trip 
through  the  cavity.  For  net  gain  P1P2GH  <  1,  the  quantity  inside  the  square 
brackets  may  be  analytically  summed,  so  that  the  output  may  be  written 


r  =  (7i72)gfi  (.s 

l-{prP2GHy 

Equation  (5)  has  the  form  of  the  Airy  formula  obtained  via  the  standard 
treatment  of  a  resonant  optical  amplifier  found  in  optics  texts  [4].  Identi¬ 
fying  pi,  P2  and  7i,72  as  the  mirror  amplitude  reflection  and  transmission 
coefficients,  respectively,  and  setting  G  =  H  =  as  the  gain  and 

phase  experienced  by  a  field  for  a  one-way  pass  through  the  cavity  of  length 
I  and  index  of  refraction  n,  with  c  the  speed  of  light  in  vacuum,  the  familiar 
transmission  characteristic  for  a  Fabry-Perot  resonant  amplifier  is  obtained: 


(6) 


1  -  PlP2G2e-^2na;7c- 

.  The  power  transmission  function  is  the  magnitude-squared  of  this  expression: 


Pputi^^)  _  _ (7i72<^o)^ _ 

Pin{uj)  ~  1  +  {piP2Giy  -  2P1P2GI  cos{2ujnl/c)  ■ 

The  resonant  amplifier  power  transmission  function  is  plotted  in  Figure 


2  for  facet  reflectivities  pi  =  p2  =  0.9  and  Go  =  1.  The  system  is  resonant 
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at  frequencies  corresponding  to  integer  multiples  m  =  1,  2,  •  •  •  of  the  inverse 
round-trip  time  through  the  cavity  r  =  2nZ/c: 


2'nm 


— 


(8) 


Equation  (5)  is  the  general  form  of  the  input-to-output  transfer  function 
for  a  delay-type  feedback  system  as  a  function  of  frequency.  This  expres¬ 
sion  is  valid  for  a  system  with  no  multiplicative  noise,  acting  on  sinusoidal 
input  fields.  As  the  feedback  fraction  is  increased,  the  widths  of  the  trans¬ 
mission  peaks  given  by  equation  (7)  decrease.  The  half-power  amplification 
bandwidth  of  a  resonance  peak  is  [4] 


Plp2Gi  ITT 

So,  as  the  net  round-trip  gain  approaches  unity,  the  bandwidth  of  the  reso¬ 
nant  amplifier  approahes  zero. 


2.2  Linear  Feedback  System  Analysis  with  Multi¬ 
plicative  Noise 

In  this  section,  a  general  formalism  will  be  developed  for  deriving  the  sys¬ 
tem  output  field  C  for  the  resonant  amplifier  of  Figure  1  when  the  forward 
and  reverse  gains  are  perturbed  by  multiplicative  noise.  The  multiplicative 
noise  may  be  fundamental  or  technical  in  origin,  as  discussed  previously.  In 
the  general  case  that  G  and  H  have  non- negligible  gain  and  phase  fluctua¬ 
tions,  the  closed-form  summation  leading  to  equation  (5)  can  not  be  applied. 
Instead,  the  partial  waves  must  be  explicitly  summed  to  obtain  the  output 


at  each  desired  time  instant  t,  and  the  amplitude  and  phase  of  the  output 
field  at  the  optical  frequency  Wm,  will  fluctuate  in  time.  This  analysis  will 
compute  the  slowly-rvarying  complex  envelope  of  the  output  field  at  a  cavity 
resonance  Um,  or 

It  is  desired  to  analyze  the  fluctuations  of  the  output  amphtude  and  phase 
of  C(urm,  t)  on  time  scales  that  are  long  compared  to  the  period  of  the  optical 
frequency  Um-  As  will  be  described,  this  can  be  conveniently  done  by  moving 
to  a  discrete-time  picture  of  the  system,  breaking  up  time  into  discrete  incre¬ 
ments  equal  to  the  time  delay  for  one  round-trip  through  the  cavity,  r.  After 
the  completion  of  this  work,  a  similar  approach  was  discovered  in  Reference 
[7]  in  which  Z-domain  techniques  from  digital  systems  theory  were  used  to 
obtain  the  transmission  and  reflection  characteristics  of  Fabry-Perot  etalons 
versus  the  optical  input  frequency.  However,  in  that  work,  no  time-variation 
of  the  feedback  coefficients  was  considered.  Instead,  the  Z-domain  analysis 
was  used  to  simphfy  the  difficult  problem  of  computing  the  transmission  and 
reflection  characteristics  of  complicated  multi-layer  structures.  In  the  present 
work,  a  similar  discrete-time  approach  is  used  on  a  simple  two-mirror  Fabry- 
Perot  structure,  but  extended  to  the  case  of  time-varying  gain  coefficients  to 
compute  the  noise  properties  of  the  output  field. 

For  the  special  case  of  a  feedback  system  operated  near  resonance,  such 
as  a  Fabry-Perot  optical  cavity,  the  throughput  is  only  appreciable  in  narrow 
bands  of  frequencies  near  the  cavity  resonances  which  occur  at  firequencies 
u>jri,  as  seen  in  the  previous  section.  All  other  frequencies  are  heavily  attenu¬ 
ated  .  In  many  practical  appUcations,  such  as  the 

laser,  the  fluctuations  of  the  cavity  transfer  function  near  a  resonance  fre- 
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quency  axe  of  primary  interest,  since  they  determine  the  long-term  stability 
of  the  laser  frequency.  Also,  the  system  operates  at  a  single  cavity  resonance, 
so  a  full  frequency  domain  analysis  is  not  really  necessary,  and  it  is  sufficient 
to  analyze  the  effect  of  the  gain  coefficient  fluctuations  on  the  behavior  of 
the  cavity  near  the  resonance  frequency  of  interest. 

Now,  the  resonant  amplifier  of  Figure  1  will  be  analyzed  near  a  single 
resonance  frequency  when  the  gain  coefficients  G  and  H  are  pertmbed  by 
noise.  To  keep  the  analysis  simple  and  the  results  transparent,  initially  it 
will  be  assumed  that  the  input  state  R  is  a  pure  sinusoid  coincident  with  a 
resonance  frequency,  Um,  with  constant  amplitude  Ar,  as  defined  in  equation 
(3).  Later,  the  more  realistic  case  of  an  input  field  with  additive  noise  will  be 
treated.  The  input  and  output  couplers  1  and  2  are  assumed  to  be  identical, 
so  that  Pi  =  P2  =  P  and  7i  =  72  =  7-  It  is  assumed  that  the  cavity 
length  I  is  time-invariant,  as  is  the  nominal  refractive  index  n.  All  phase 
or  gain  fluctuations  are  described  by  complex  zero-mean  random  variables 
6g{t)  and  6h{t).  The  time- varying  forward  and  reverse  gain  coefficients  due 
to  multiplicative  noise  are  then  modeled  as 

t)  =  Go  G{t).  (10) 

H{um,  t)  =  Ho  =  Ho  H{t).  (11) 

where  Go,  Ho  represent  the  static  forward  and  reverse  loss  or  gain, 
is  the  static  one-way  phase  shift  added  to  a  signal  at  frequency  Um  passing 
through  the  cavity.  At  a  loop  resonance  Wm,  the  static  phase  shift  will  satisfy 
the  condition  nuJmllc  =  27rm  radians,  with  m  an  integer.  6g(t)  and  8h{t)  are 
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complex  random  variables  representing  the  time-varying  phase  and  gain  at 
frequency  oJm-  These  may  be  written 


6g{t)  =  6g'{t)  -h  i5g"{t) 

(12) 

Shit)  =  6h!{t)  +  i6h"it) 

(13) 

where  the  real  parts  Sg'(t)  and  6h'{t)  represent  phase  fluctuations,  and  the 
imaginary  parts  6g"{t)  and  8h"{t)  represent  gain  fluctuations.  The  statistical 
description  and  any  correlation  of  the  real  and  imaginary  parts  will  depend 
in  detail  on  the  physics  of  the  gain  medium. 

The  transformation  from  a  continuous-time  representation  to  the  discrete¬ 
time  representation  is  accomplished  by  considering  the  static  phase  shifts 
g-inoimi/c  as  delay  elements  of  duration  r/2,  as  illustrated  in  Figure  3,  and  rep)- 
resenting  the  continuous-time  gain  variations  of  G{t)  and  H{t)  as  a  discrete¬ 
time  random  sequence  of  complex  values  Gt  and  Ht  with  dc  gain  Go  and 
Ho-. 


Gt  =  Go  =  Go  Gt.  (14) 

Ht  =  Ho  =  Ho  Hf  (15) 

separated  by  time-intervals  r.  (Throughout  the  rest  of  the  analysis,  discrete¬ 
time  variables  will  be  denoted  with  a  subscripted  time  index,  e.g.,  Xt,  to 
distinguish  them  from  continuous-time  quantities;  X{t)).  In  this  discrete¬ 
time  picture,  C{uJm,t)  is  assumed  to  be  a  pure  sinusoid  at  loop  resonance 
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frequency  a;^,  with  slowly- varying  amplitude  and  phase  fluctuations,  which 
will  be  referred  to  as  C*. 

Now,  we  apply  this  discrete-time  definition  of  the  gain  to  the  linear  anal¬ 
ysis  of  the  previous  section  to  obtain  an  expression  for  the  output  field  Cf 
At  time  increment  t,  the  forward  gain  coefficient  is  Gt,  the  reverse  gain  is  Ht, 
and  the  input  state  is  Rt-  For  notational  clarity,  it  is  understood  throughout 
that  time  index  t-  j  refers  to  time  t  -  jr,  since  time  has  been  quantized  in 
units  of  the  round-trip  time  in  the  cavity  r.  The  output  state  is  then  written 

Ct  =  'y^GtRt  +  'y^p^GtHt-iGt-iRt-i 

p‘^GtHt-iGt-iHt-2Gt-2Rt-2  H - 

=  'y'^Gt[Rt  +  Ht-\Gt-\Rt-i 

+p^Ht-iGt-iHt-2Gt-2Rt-2  +  •  •  •]  (16) 

It  is  assumed  that  the  dc  forward  and  reverse  gains  and  static  phase  are 
identical,  so  that  Go  =  Ho-  This  assumption  is  equivalent  to  requiring  the 
loss  (or  gain)  and  time-of-flight  to  be  equal  for  both  directions  of  propagation 
in  the  optical  cavity.  This  is  a  good  assumption  for  a  Fabry-Perot  cavity, 
since  the  light  passes  through  the  same  physical  medium  in  both  directions, 
but  may  not  be  valid  for  other  cavity  geometries.  It  is  further  assumed  that 
the  gain  fluctuation  rate  is  slow  compared  to  the  time-of-flight  r  for  a  round- 
trip  in  the  cavity.  Then,  it  is  valid  to  assume  the  forward  and  reverse  gain 
fluctuations  are  equal  over  any  round-trip  time  interval  r,  so  the  round-trip 
gain  may  be  written  GtHt  =  H^.  Then,  defining  the  time-varying  net  round- 
trip  gain  to  be  p'^H^  =  Kt,  the  output  field  at  time  t  is  given  by  the  infinite 


series: 


Ct  =  Ht[Rt  +  Kf-iRt-i  +  Kt-iKt-2Rt-2  H - ]•  (17) 

To  analyze  the  characteristics  of  the  gain  coefficient  fluctuations  on  the 
output  fleld,  first  the  case  of  an  ideal  noiseless  input  field  will  be  examined. 
The  input  state  is  taken  to  be  a  constant  Rt  =  R  for  all  times  t,  so  it  may 
be  factored  out  of  equation  (17).  Then,  the  output  field  Ct  is  expressible  as 
an  infinite  sum  of  partial  waves  in  the  cavity  at  any  given  time  instant  t: 


Ct  =  7^HtR[l  +  Kt-i 


+  {Kt-iKt-2)  +  {Kt-\Kt-2Kt-z)  H - ]. 


(18) 


This  infinite  series  will  be  truncated  to  a  finite  number  of  M  terms,  as 
follows.  If  the  gain  fluctuations  are  small,  then  the  amplitudes  of  the  succes¬ 
sive  terms  are  determined  by  increasing  powers  of  the  static  net  round-trip 
gain  coefficient,  so  that  the  j  -  1st  term  is  proportional  to  Ki.  When  the 
mean  value  of  the  net  round-trip  gain  is  less  than  unity,  Ko  =  <  1,  the 

infinite  sum  of  Ki  is  equal  to  l/{l-Ko),  as  was  used  previously  to  derive  the 
Airy  formula.  If  a  finite  number  of  terms  are  used  instead,  the  summation 
of  Ki  from  j  =  0  to  j  =  M  is  given  by 


M 


j=0 


1  -  K^+^ 
l-Ko 


(19) 


The  value  of  M  is  determined  by  requiring  the  difference  between  the  infinite 
sum  and  the  finite  sum  of  M  -i- 1  terms  to  be  less  than  a  small  fraction  e  of 
the  value  of  the  full  infinite  sum.  This  condition  is  written 
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(i^)  -  (t^)  <  ^  (t^ 

Then,  solving  for  the  required  number  of  terms  M  yields 


(20) 


"  >  5Si  -  >’» 

We  shall  require  e  <  0.01,  or  the  difference  between  the  finite  sum  and  the 
infinite  sum  of  less  than  1  percent.  For  example,  if  Ko  =  0.9,  then  from 
equation  (21)  we  obtain  M  >  43.  The  time  interval  Mr  is  proportional  to 
the  memory  of  the  feedback  system,  and  increases  as  the  net  round-trip  gain 
K  approaches  unity.  However,  it  will  be  seen  later  that  the  value  of  M  just 
derived  is  only  valid  in  the  case  of  small  gain  fluctuations. 

The  output  field  at  any  time  t  is  a  random  phasor  sum  of  the  M  -f  1 
partial  waves.  This  expression  may  be  cast  in  a  more  illuminating  form  by 
substituting  the  full  form  of  the  discrete  gain  coefficients  for  the  Kt  terms 
given  by  equation  (15).  Doing  this  and  simplifying  yields: 


Ct  =  +  p'Hl  -t- 

_|_  .  .  .  -j_  g— - (22) 

Each  successive  term  of  this  expression  represents  an  individual  partial  wave 
that  has  travelled  an  increasing  number  of  times  through  the  cavity,  and  has 
thus  experienced  the  gain  and  phase  fluctuations  6ht  over  an  increasing  time 
span.  For  the  case  of  no  gain  or  phase  fluctuations,  (i.e.,  6ht  =  0  for  all  t,) 
this  expression  may  be  summed  analytically,  and  collapses  to  the  well-known 
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Airy  formula  for  the  Fabry- Perot  transmission  derived  previously.  But  in  the 
present  form,  each  partial  wave  has  a  diflFerent  random  gain  and  phase,  due 
to  the  multiplicative  complex  fluctuations  8ht. 

Prom  equation  (22)  it  is  seen  that  the  phases  of  the  successive  partial 
waves  are  cumulative  sums  of  the  multiplicative  noise  process  8ht  over  in¬ 
creasing  times.  If  8ht  is  a  white  Gaussian  complex  noise  process  with  zero 
mean  and  variance  then  the  complex  phases  of  the  successive  terms 
represent  a  random  walk  process.  The  variance  of  a  random-walk  process 
grows  linearly  with  time,  so  the  complex  phase  of  the  Mth  partial  wave, 
2{8ht-i  +  8ht-2  H - +  8ht-M),  wiU  have  variance  alf  =  2Ma\. 

There  are  two  regimes  which  we  may  identify  with  respect  to  the  phase¬ 
spreading  of  the  partial  waves  in  the  output  field  summation  of  equation  (22), 
which  we  shall  denote  as  the  small-angle  and  large-angle  regimes.  The  small- 
angle  regime  corresponds  to  relatively  weak  multiplicative  noise  and/or  low 
net  round-trip  gain,  such  that  the  root-mean-square  complex  phase  of  the  last 
(i.e.,  (M-|-  l)st)  partial  wave  is  less  than  0.1.  In  this  case,  the  exponentials  in 
equation  (22)  may  be  linearized  by  the  small-angle  approximation,  and  the 
summation  for  Ct  may  be  obtained  analytically,  as  will  be  shown. 

The  large-angle  regime  corresponds  to  the  case  when  the  the  rms  com¬ 
plex  phase  variance  of  the  older  partial  waves  exceeds  0.1.  In  this  case,  the 
small-angle  approximation  does  not  obtain,  and  the  full  phasor  expression 
of  equation  (22)  must  be  used  to  compute  the  output  field  Ct  at  each  time 
increment  t.  For  a  large  multiplicative  noise  variance  cr^,  and/or  high  net 
round-trip  gain  A'o  — >  1,  the  random- walking  phase  fluctuations  of  the  later 
partial  waves  may  become  large  enough  to  cause  destructive  interference  with 
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the  earlier  partial  waves,  causing  large  fluctuations  of  the  output  fleld  am¬ 
plitude  and  phase.  Also,  the  random-walking  gain  fluctuations  of  the  later 
partial  waves  will  eventually  cause  significant  deviations  from  the  exponen¬ 
tial  decay  of  the  partial  waves  given  by  Kl.  This  causes  the  value  of  M  to 
increase,  which  leads  to  enhanced  low-frequency  fluctuations  of  the  output 
field. 

First,  the  small-angle  regime  will  be  analyzed.  In  this  regime,  the  phase 
fluctuations  of  the  output  field  due  to  multiplicative  phase  noise  are  en¬ 
hanced  as  the  net  round-trip  gain  is  increased.  Then,  the  large-angle  regime 
will  be  examined  by  computing  the  full  phasor  expression  of  equation  (22) 
for  constant  multiplicative  phase  noise  variance  while  progressively  increas¬ 
ing  the  net  round-trip  gain.  The  transition  from  the  small-angle  to  large 
angle  regimes  is  marked  by  the  onset  of  rapid,  large  fluctuations  of  the  out¬ 
put  field  amplitude  and  phase,  and  a  decrease  in  the  average  output  power. 
The  results  obtained  will  be  used  to  calculate  the  linewidth  of  the  resonant 
amplifier  in  each  regime.  In  the  small-angle  regime,  it  is  found  that  the  mul¬ 
tiplicative  phase  noise  leads  to  an  output-power-independent  contribution  to 
the  linewidth  of  the  output  field  power  spectrum.  In  the  large-angle  regime, 
the  linewidth  increases  with  increasing  net  gain.  Both  of  these  results  are 
counter  to  the  decreasing  linewidth  with  increasing  net  gain  expected  from 
equation  (9)  for  the  standard  analysis  without  multiplicative  noise. 
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2.3  Small-angle  multiplicative  noise  regime:  linear 
approximation 

In  this  section,  the  case  of  relatively  weak  multiplicative  noise  will  be  treated 
using  a  linear  approximation.  The  last  term  of  equation  (22)  represents  the 
“oldest”  partial  wave  that  has  travelled  M  +  1  round-trips  in  the  cavity.  The 
complex  phase  of  this  wave  is  an  M  -1-  1-step  random  walk.  The  small-angle 
regime  is  defined  as  the  case  when  the  standard  deviation  of  the  complex 
phase  for  the  M  -H  1st  partial  wave  is  less  than  0.1,  or 

Vmcxh  <  0.1  (23) 

where  M  is  defined  by  equation  (21).  Then,  the  small-angle  approximation 
for  the  exponential  (e®  l-t-  x)  can  be  used  to  linearize  the  expression  for 
the  fluctuating  output  field  so  that  equation  (22)  may  be  rewritten: 

Ct  «  7"i?o(l  +  *(<5/it))ii:[l  +  P^Hl  (1  -f-  2z((5/it_i))  -H  (1  +  2i{8ht-i  +  <5/it-2)) 
+  •  •  •  +  p^MHl^  (1  +  2i{8ht-i  +  6ht-2  +  •  ■  •  +  Sht-M))]-  (24) 

Substituting  Ko  =  p^H^  into  this  expression  yields: 

Ct  =  +  i{6ht))R[l  +  Ko{\  +  2i{6ht-x))  +  K]  (1  -i-  2i{8ht-i  +  8ht.2)) 

(1  2i{8ht-\  -h  8ht-2  +  •  •  ■  +  8ht-M))]-  (25) 

Expanding  this  expression  will  yield  many  terms  of  Oifih?),  like  {8h^8ht-k), 
which  are  small  compared  to  0{8h),  and  so  may  be  discarded,  leaving 
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G  =  ^^HoR[{l  +  Ko  +  K^o+---  +  K^)  +  i^  +  Ko  +  Kl  +  ---  +  K^){i6ht) 

+{Ko  +  KI  +  ---  +  K^){2i6ht-.,)  +  {K^,  +  +  . . .  +  K^){2i6ht^,) 

+  -.-  +  K^{2i6ht^M)].  (26) 

The  term  proportional  to  6ht  is  due  to  the  phase  fluctuation  encountered 
in  the  first  “half  round-trip”  through  the  cavity.  All  other  terms  represent 
“full  round-trips”  through  the  cavity,  hence  the  factors  of  2  appear  in  terms 
Sht-i,6ht-‘2,  •  •  •  etc. 

The  summations  of  powers  of  Ko  may  be  performed  analytically  using 
equation  (19).  The  coefficients  in  equation  (26)  are  then  designated  as  ak, 
where  the  coefficient  of  the  first  term  is  oo,  defined  as 

(M+l)-l  1  _  t^M+1 

“»=  E  Ki=  i_y-" 

j=0  ^ 

and  subsequent  terms  a^,  A:  >  1  may  be  written 


j=0 

1  -  1  - 

I -Ko  I -Ko 

K^o  - 

l-Ko  ■ 

As  M  — >  oo,  Uik  converges  to 


1-K„ 


(28) 


(29) 
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Now,  substituting  the  coefficients  into  equation  (26),  for  the  present 
special  case  of  a  noiseless  input  R,  the  output  field  at  any  time  instant  t'is 
seen  to  be  a  weighted  summation  of  the  multiplicative  noise  samples  ht  over 
the  time  interval  t  —  M  to  t: 


Ct  —  'y^  HoR[ao  +  aoiSht  +  a,\i26ht-x  +  ■  •  •  + 


=  I'^HoR 


M 


ao(l  +  iSht)  +.  ^  flfc  i25ht-k 


fc=i 


(30) 


Finally,  the  factor  1/(1  —  Ko)  may  be  factored  out  of  this  expression  to  yield 


a  =  [l  +  iSht  +  iKo26ht.x  +  iKl26ht.2  +  •  •  •  +  iK^26ht-M\  ■ 

i  £\q  ^ 

(31) 

This  expression  represents  the  transformation  of  the  multiplicative  phase 
and  gain  fluctuations  to  the  output  field  at  every  time  instant  t.  The  first 
term  is  the  static  gain  of  the  resonant  amplifier  as  given  from  the  standard 
analysis  leading  to  equation  (5).  Recall  from  equation  (13)  that  6ht  is  a 
complex- valued  random  series  whose  real  and  imaginary  parts  represent  the 
time-varying  gain  and  phase,  respectively,  of  the  amplifying  medium.  Phys¬ 
ically,  the  real-valued  coefficients  Uk  determine  the  contribution  of  the  mul¬ 
tiplicative  phase  and  gain  fluctuation  at  the  previous  time  increment  t  ^  k 
to  the  output  at  time  increment  t.  It  is  apparent  from  equation  (31)  that 
for  \Ko\  <  1,  the  a*,  magnitudes  decrease  exponentially  with  increasing  k. 
Thus,  noise  events  further  removed  in  time  from  the  time  increment  t  have 
a  diminishing  effect  on  the  output  state.  In  this  small-angle  approximation. 
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noise  events  occurring  at  times  prior  to  t  —  M  have  no  effect  on  the  output, 
which  is  just  an  analytic  statement  of  the  finite  memory  time  of  the  system. 

In  the  analysis  thus  far,  fluctuations  of  the  input  field  R  were  not  consid¬ 
ered.  In  general,  the  input  R  =  Rt  may  also  have  a  time-varying  amplitude 
and  phase.  To  treat  this  possibility,  the  input  state  is  now  defined  in  the 
continuous-time  picture  as  a  sinusoid  with  mean  frequency  coincident  with 
a  cavity  resonance  Um,  perturbed  by  amplitude  noise  6A{t)  and  phase  noise 
64>{t) 

R{t)  =  {Ar  +  6A{t))  (32) 

Assuming  S(t)t  <0.1  radian,  the  small-angle  approximation  can  again  be 
employed,  and  the  time-varying  input  field  R{t)  may  be  manipulated  into 
the  discrete-time  form: 

Ri  {Ar  +  6At)  (1  -h 

=  {Ar  -1-  iAR6(j)t  -h  6 At  +  i5At6(j)t)e^‘^”'^ 

Ar  J 

=  ARe^-'^l  +  6Rt) 

Rt  =  Ro{l  +  6Rt),  (33) 

where  the  term  of  0{6A){6^)  was  neglected,  since  it  is  small  compared  to 
the  terms  of  0{6A)  and  0{5(l)),  and  the  noise  sources  6 A  and  Scf)  are  assumed 
to  be  uncorrelated. 

Now,  this  expression  for  the  fluctuating  input  field  is  inserted  into  the  full 
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form  of  equation  (22)  to  obtain  the  output  at  every  time  increment  t  due  to 
both  input  field  fluctuations  and  multiplicative  gain  fluctuations: 


p^GtHt-iGt-iHt-2Gt-2Rt-2  +  •  •  ■ 

=  'y^GtlRt  +  Ht-iGt-iRt-i 

-\-p‘^Ht-iGt-iHt-2Gt-2Rt-2  H - ] 

=  'y'^Gt[Rt-\-  Kt-iRt-i-^  Kt-\K.t-2Rt-2-\-  •  •  •]■  (34) 

Using  manipulations  similar  to  those  used  to  derive  equation  (31)  (the  de¬ 
tailed  algebra  is  relegated  to  an  Appendix)  a  simplified  expression  for  the 
complex  output  field  amplitude  is  obtained: 

Ct  =  [1  +  iSht  +  Koi25ht-iKli26ht-2  +  •  •  •  +  K^i26ht-iM-i) 

i  -^O 

+(1  -  Ko)  {5Rt  +  Ko6Rt-i  +  K^JRt-2  +  K^Rt-3  +  •  •  •  +  /(Tf  (5i?t_^^) 

The  powers  of  Ko  multiplying  the  6R  terms  may  be  defined  as  a  set  of 
coefficients  bk'- 


h  =  Kl  (36) 

The  full  expression  for  the  output  field  can  now  be  written  more  compactly  as 
summations  of  the  two  noise  processes  6ht  and  6Rt  in  terms  of  the  coefficients 
Ufe  and  bk'. 
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Ct  =  j^HoR  ao(l  +  Sht)  +  aki2Sht-.i}\  +  (Y^bkSRt-i^  .  (37) 

This  expression  represents  a  linear  approximation  to  the  full  form  of  equation 
(22)  for  the  slowly-varying  complex  envelope  of  the  output  field  Ct  including 
both  input  field  fluctuations  and  multiplicative  phase  and  gain  fluctuations. 
This  linear  approximation  is  valid  in  the  small-angle  regime  as  defined  by 
equation  (23). 

Now,  the  power  spectrum  of  the  output  field  in  the  small-angle  regime 
will  be  calculated.  Equation  (37)  has  the  form  of  a  finite-impulse-response 
filter  operating  on  “input”  processes  6ht  and  SRf.  For  linear  systems  in 
general,  the  Fourier  transform  of  the  time  impulse-response  yields  the  output 
frequency  response  of  the  system  [8].  For  an  input  field  noise  source  modeled 
by  a  Gaussian-distributed  random  process,  the  power  spectrum  of  the  output 
field  is  the  product  of  the  input  noise  power  spectrum  and  the  magnitude- 
squared  of  the  system  frequency  response  [?].  So,  if  the  impulse  responses, 
net  round-trip  gain  Kg,  and  the  analytical  form  of  the  power  spectra  of  the 
multiplicative  fluctuations  Sh{f)  and  the  input  field  fluctuations  5'/{(/)  are 
known,  then  the  output  power  spectrum  Sc{f)  may  be  obtained  analytically. 

The  impulse-responses  of  the  resonant  amplifier  output  to  a  delta-function 
of  multiplicative  noise  6ht  or  input  noise  6Rt  are  defined  to  be  6Ch,t  and  SCR^t, 
respectively.  The  multiplicative  impulse-response  6Ch,t  is  the  decaying  out¬ 
put  field  response  given  by  equation  (37)  due  to  a  delta-function  multiplica¬ 
tive  impulse  in  6ht  applied  to  the  system  at  time  t  =  0: 


5Ch,t  = 


'fHoR  2at  Ut 


(2f^\ 

W-Ko)) 


2j^HoR 
(1  -  Ko) 


(38) 


where  the  decay  constant,  and  Ut  is  the  unit  step  function,  with  Ut  =  0  for 
t  <  0,  and  Ut  =  1  for  t  >  0,  and  it  is  assumed  that  t  takes  only  integer 
values  (i.e.,  time  is  normalized  to  increments  of  r).  The  first  half-round-trip 
(/c  =  0  term)  was  taken  to  be  equal  to  a  full  round-trip.  For  M  >  100 
corresponding  to  Kq  >  0.95,  this  introduces  a  small  error  of  less  than  a  1 
percent  in  the  output  sum.  The  response  to  an  input  field  impulse  in  6Rt 
is  obtained  similarly  as  the  response  of  equation  (37)  to  a  delta-function 
impulse  in  8Rt  applied  to  the  system  at  time  t  =  0: 


8CR^t  =  l^HoR  bt  Ut 
=  ^'^HoRKlut 

=  'i^EoR  Uf  (39) 

The  Fourier  transforms  of  the  impulse  responses  8Ch,t  and  bCR^t  yield  the 
corresponding  frequency  responses  F/i(/)  and  Ffi(/),  where  /  is  the  Fourier 
frequency.  Both  impulse  responses  represent  decaying  exponentials,  leading 
to  Lorentzian  frequency  responses  given  by 


FM  = 


(  2'i^EoR  \ 

Ui-^o)ry 


1 


1  -i- 


(40) 
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and 


FrU)  = 


(^) 


(41) 


where  Uc  =  T/t  is  the  rolloff  frequency  of  the  Lorentzian  response. 

The  power  spectrum  of  the  output  field  is  comprised  of  three  terms:  a 
delta-function  term  due  to  the  sinusoidal  input  R,  and  terms  due  to  the  mul-^ 
tiplicative  and  input  fiuctuations.  In  the  absence  of  multiplicative  noise,  and 
for  a  pure  sinusoidal  input  field  i?  at  a  resonance  frequency  Wm,  the  output 
power  spectrum  Sc{f)  is  a  pure  sinusoid  of  frequency  Um,  with  magnitude 
given  by  the  constant  term  in  equation  (37).  This  may  be  written  as 


Seif)  = 


where  the  delta  function  is  defined  as  6(0)  =  1  and  6{x)  =  0,x  ^  0,  and 
is  centered  on  frequency  Um-  The  Fourier  firequency  /  is  the  offset  in  Hz 
from  the  optical  frequency  This  is  equivalent  to  the  power  gain  at 
the  peak  of  the  resonant  amplifier  transmission  as  given  by  equation  (7). 
The  noise  processes  6ht  and  SRt  give  rise  to  additional  terms  in  the  output 
power  spectrum.  It  is  assumed  that  the  noise  processes  6ht  and  6Rt  are 
uncorrelated,  so  the  power  spectra  due  to  each  process  may  be  added.  The 
term  due  to  multiplicative  noise  is  the  product  of  the  multiplicative  noise 
power  spectral  density  Sh{f)  and  the  magnitude-squared  of  the  frequency 
response  Fh{f)-  Similarly,  the  term  due  to  the  input  noise  is  the  product  of 


{j^HoRaof  6{f) 

m'x" 


(42) 
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the  input  field  noise  power  spectral  density  S'i?(/)  and  the  magnitude-squared 
of  the  frequency  response  The  addition  of  equation  (42)  and  the  two 

noise  terms  yields  the  total  output  power  spectrum: 

Sc(f)  =  W)j  +  ie.(/)t  S»(/)  +  Mf)  (43) 

Substituting  for  the  frequency  responses  Fhif)  and  Fh(/)  and  subsequent 
simplification  yields 


ScU)  = 


Shif)  + 


1  +  i 


2irf 

OJc 


+  TtMt 


VI 

Consider  the  last  term  of  this  expression  due  to  the  input  noise  bR.  For 
low  gain  and/or  low  multiplicative  noise,  the  signal-to-noise  ratio  inside  the 
amplifier  bandwidth  will  be  the  same  as  the  input  signal-to-noise  ratio.  Phys¬ 
ically,  the  last  term  implies  that  the  output  field  phase  and  amplitude  exactly 
follows  the  input  field  phase  and  amplitude  fiuctuations  for  rates  less  than 
the  cavity  bandwidth  azc,  which  is  intuitively  correct.  For  fluctuation  rates 
faster  than  the  input  field  fluctuations  are  rolled  off  as  1//^,  as  expected, 
since  they  lie  outside  the  bandpass  of  the  cavity.  As  the  net  gain  Ko  is 
increased,  Uc  decreases,  and  the  integrated  noise  power  due  to  the  input  fluc¬ 
tuations  decreases,  implying  a  decreased  spectral  width  of  the  output  field. 


S«(/) 

(44) 


B.34 


The  resonant  amplifier  therefore  produces  a  narrowband  output  signal  of 
high  spectral  purity  from  a  broadband  noise  input  when  the  net  round-trip 
gain  Ko  is  high.  This  is  the  expected  behavior  of  a  resonant  amplifier  as 
obtained  from  the  standard  analysis  leading  to  equation  (7),  and  is  the  basis 
of  the  linear  model  of  the  laser  as  a  noise  driven  resonant  amplifier  [4],  which 
we  will  return  to  later. 

Now  consider  the  second  term  of  equation  (44)  .  It  is  seen  that  the  second 
term  due  to  multiplicative  noise  6ht  is  multiplied  by  the  factor  1/(1  -  KoY 
compared  to  the  last  term  due  to  the  input  field  fluctuations  6Rt.  What 
is  significant  about  this  is  that  as  the  net  gain  Ko  is  increased,  given  fixed 
levels  of  input  fluctuations  5fi(/)  and  multiplicative  fluctuations  Sh{f),  even¬ 
tually  the  output  field  power  spectrum  Sc{f)  will  become  dominated  by  the 
multiplicative  noise.  After  this  happens,  the  signal-to- noise  ratio  for  Fourier 
frequencies  less  than  Uc  will  decrease  with  further  increases  in  the  net  gain. 

Making  use  of  the  fact  that  for  TCo  ~  1,  T  =  \^n{Ko)\  ~  (1  —  Ko),  the 
rolloff  frequency  may  be  approximated  as  «  (1  —  Ko)/t.  Making  this 
substitution  in  the  second  term  of  equation  (44)  yields 


ScU)\ 


Ko^l 


V-Ko^ 


S{f)  + 


(1- 7^0)2 +  (27r/r)^ 


Shif)  + 


1 


SRifm 


In  the  limit  as  ATo  1,  the  second  term  due  to  multiplicative  noise  over¬ 
whelms  the  last  term,  so  the  output  power  spectrum  becomes 
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(46) 

In  the  limit  a.s  Kg  I,  the  integrated  multiplicative  noise  power  given  by 
the  second  term  in  equation  (45)  approaches  a  constant  value,  which  implies 
that  the  spectral  width  of  the  output  signal  also  becomes  constant.  This  is  in 
marked  contrast  to  the  prediction  of  the  standard  analysis  [4]  of  equation  (7), 
which  predicts  decreasing  linewidth  with  increasing  gain.  In  summary,  for 
low  net  gain  when  the  input  field  fluctuations  are  dominant,  the  width  of  the 
output  power  spectrum  decreases  with  increasing  gain.  But  as  the  net  gain 
Ko  1,  the  multiplicative  fluctuations  eventually  become  dominant,  and 
the  power  spectrum  of  the  resonant  amplifier  approaches  a  constant  value. 
Therefore,  a  minimum  linewidth  should  be  expected  for  resonant  amplifiers 
operated  at  high  gain  in  the  presence  of  multiplicative  noise.  To  the  extent 
that  the  multiplicative  fluctuations  arise  from  fundamental  thermodynamic 
processes  [12],  this  minimum  linewidth  will  be  a  fundamental  limit. 

2.4  Large-angle  multiplicative  noise:  numerical  anal¬ 
ysis 

In  this  section,  the  full  phasor  form  of  the  output  field  given  by  equation 
(22)  is  analyzed  in  the  large-angle  regime,  by  numerical  simulation  using 
computer-generated  random  time  series  to  represent  the  multiplicative  noise. 
This  approach  is  taken  to  investigate  the  behavior  of  the  output  field  ampli¬ 
tude  and  phase  when  the  variance  of  the  complex  phase  of  the  older  partial 
waves  exceeds  0.1  radians,  so  the  small-angle  approximation  which  permits 
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the  analytic  derivation  of  equation  (23)  no  longer  obtains.  The  output  field 
time  series  is  computed  for  a  noiseless  input  field  R  of  unit  amplitude  at  fixed 
cavity  resonance  frequency  Um,  for  various  values  of  the  net  round-trip  gain. 
The  calculations  are  performed  for  a  300//m-long  Fabry-Perot  semiconduc¬ 
tor  laser  cavity,  at  1.3/im  wavelength,  with  equal  facet  power  reflectivities 
=  0.32. 

The  real  and  imaginary  parts  of  6ht  are  assumed  to  be  correlated,  as  is 
the  case  for  phase  and  gain  fluctuations  arising  from  electron-number  fluc¬ 
tuations  in  a  semiconductor  laser  gain  medium  [10].  The  real  and  imaginary 
parts  of  Sht  are  written  in  terms  of  a  common  noise  source  6nt  representing 
fluctuations  of  the  refractive  index: 

6h{t)  =  6h'{t)  -I-  i6h"{t) 

The  constant  P  relates  the  changes  in  the  real  and  imaginary  parts  of  the 
refractive  index  due  to  electron  density  fluctuations,  and  is  a  measure  of 
the  amount  of  amplitude-phase  coupling  for  waves  in  the  gain  medium.  It 
was  shown  [2],  [3]  that  P  is  responsible  for  broadening  of  the  semiconductor 
laser  linewidth  above  the  Schawlow-Townes  prediction,  and  so  P  is  typically 
referred  to  as  the  linewidth  enhancement  factor.  The  refractive  index  in  a 
semiconductor  laser  gain  medium  is  typically  modeled  as  a  linear  function  of 
the  carrier  density  [?].  Thus,  changes  in  the  electron  density  due  to  injection 
current  noise,  thermal  fluctuations,  or  electron-hole  recombination  processes 
cause  correlated  changes  in  the  gain  and  phase  of  the  amplifying  medium  via 


Computer  simulations  were  run  for  increasing  net  round-trip  gains  Ko  = 
0.9,  0.95,  0.99,  0.995,  and  0.999,  which  represents  three  orders  of  magnitude 
change  in  the  net  gain.  For  each  value  of  /f©,  an  iV-element  complex  vector 
of  the  fluctuating  output  field  Ct  is  produced.  For  these  net  gain  values,  the 
total  numbers  Af  of  partial-waves  used  in  the  summations  were,  respectively, 
M  —  50,  100,  500,  1000,  and  5000.  The  complex  multiplicative  noise  process 
6ht  is  simulated  by  a  computer-generated  Gaussian-distributed  zero-mean 
random  time  series  with  standard  deviation  of  Oh>  =  0.1.  This  implies  that 
the  system  is  in  the  large  angle  regime  for  all  of  the  net  gain  values  tested. 
The  computer  produces  a  series  of  random  values  from  a  Gaussian  distribu¬ 
tion  such  that  the  rms  value  of  an  infinite  number  of  such  values  would  be 

ffv. 

The  computed  results  for  net  round-trip  gain  of  Ko  =  0.9,  0.95,  0.99, 
0.995,  smd  0.999  are  now  discussed.  To  see  the  qualitative  effect  of  the 
increasing  net  gain  on  the  output,  the  real  and  imaginary  parts  of  the  output 
electric  field  Ct  are  displayed  parametrically  in  the  complex  plane  in  Fig  2 

-  6,  for  increasing  values  of  Ko-  In  these  plots,  each  of  the  N  dots  on  the 
complex  plane  represents  the  computed  position  of  the  tip  of  the  output  field 
vector  Ct.  Figure  2  illustrates  the  case  of  net  round-trip  gain  Ko  =  0.9.  The 
tip  of  the  electric  field  vector  describes  an  arc  with  maximum  phase  excursion 
of  approximately  +/-  0.3  radians,  but  with  relatively  constant  amplitude.  As 
the  gain  is  increased  to /fo  —  9.95, 0.99, 0.995,  and  0.999,  illustrated  in  Fig  2 

-  6,  the  maximum  phase  excursion  increases  and  the  amplitude  exhibits 
progressively  larger  fluctuations  which  appear  as  a  swirling  pattern  in  the 
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complex  plane.  If  there  were  no  multiplicative  noise,  all  of  the  dots  would 
fall  on  one  point  on  the  real  axis,  determined  by  the  static  gain  1/(1  —  Ko). 
The  small-angle  regime  treated  previously  represents  small  excursions  about 
this  point. 

As  the  output  field  vector  swirls  around  the  complex  plane,  Figure  4e, 
the  magnitude  and  phase  change  drastically.  The  phase  and  instantaneous 
intensity  (magnitude-squared  of  the  electric  field)  versus  time  for  each  value 
of  net  gain  are  plotted  in  Figures  7  -  16 .  As  the  net  gain  increases,  the 
period  and  amplitude  of  both  the  intensity  and  phase  fluctuations  increase. 
At  the  lower  gain  values,  the  mean  value  of  the  intensity  fluctuations  is 
close  to  the  value  expected  without  multiphcative  noise.  However,  as  the 
gain  is  increased,  the  amplitude  of  the  intensity  fluctuations  increase  greatly. 
The  intensity  becomes  more  deeply  modulated,  and  the  mean  value  of  the 
intensity  moves  to  lower  values.  This  behavior  is  more  apparent  from  the 
probability  distribution  of  the  intensity,  discussed  below.  As  the  net  gain  is 
increased,  the  phase  fluctuation  amplitude  also  increases  ('Fig  12-16  )  until 
eventually,  the  phase  “wraps  around”  more  than  tt  radians. 

The  probabiUty  distribution  of  the  power  versus  increasing  net  gain  is 
plotted  in  Fig  17  -  2lAt  Ko  =  0.9,  the  probability  distribution  has  devel¬ 
oped  a  “tail”  stretching  to  lower  values  of  intensity.  As  the  net  gain  increases, 
this  tail  becomes  more  pronoimced,  and  the  mean  value  of  the  probability 
distribution  decreases.  At  high  gain  of  Ko  =  0.999,  the  intensity  probability 
distribution  is  peaked  near  zero,  nearly  resembling  the  probability  distribu¬ 
tion  of  a  thermal  emission  source,  or  a  below-threshold  laser.  The  shift  of  the 
mean  intensity  to  lower  values  is  in  qualitative  agreement  with  a  theoretical 
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Fokker-Planck  analysis  and  experimental  results  obtained  for  a  helium-neon 
laser  with  multiplicative  loss  noise  [15].  This  shift  of  the  probability  distribu¬ 
tion  to  lower  values  with  increasing  gain  is  contrary  to  the  behavior  expected 
for  the  resonant  amplifier  without  multipUcative  noise.  In  the  context  of  laser 
operation,  this  result  implies  that  the  mean  value  of  the  intensity  probabiUty 
distribution  will  initially  grow  as  the  gain  is  first  increased,  but  eventually 
will  reach  a  maximum  value,  and  then  begin  to  decrease. 

InFig  22  -  26the  power  spectrum  of  the  output  field  versus  increasing 
net  gain  is  plotted.  The  noise  level  progressively  increases  as  the  net  gain 
is  increased.  This  is  not  unexpected,  given  the  increasingly  large  intensity 
and  phase  fluctuations  of  the  output  field  as  the  net  gain  is  increased.  In 
the  absence  of  multipUcative  noise,  the  output  power  spectrum  would  be  a 
delta-function,  since  no  input  noise  source  is  included  in  the  simulations. 

It  is  noted  that  in  all  of  these  results,  the  variance  of  the  multipUcative 
noise  time-series  was  kept  constant,  and  only  the  net  gain  was  changed.  While 
such  a  large  value  of  multipUcative  noise  may  not  typically  encountered  in 
practice,  these  simulation  results  serve  to  iUustrate  the  effects  of  increasing 
net  gain  on  the  output  field.  By  using  a  relatively  large  value  of  multipUcative 
noise,  it  was  possible  to  study  the  qualitative  behavior  of  the  output  versus 
net  gain  using  a  manageable  number  of  partial  waves.  In  a  semiconductor 
laser,  net  gains  oi  Ko  ^  I  -  10"®  are  typical,  requiring  approximately  one 
million  partial  waves  per  time  increment  to  simulate.  To  obtain  a  reasonable 
estimate  of  the  output  power  spectrum,  at  least  1024  time  increments  need 
to  be  computed,  which  becomes  computationally  prohibitive. 
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2.5  Detailed  description  of  simulation  algorithm 

For  each  value  of  Ko,  the  details  of  the  simulation  algorithm  are  as  follows: 
First,  a  real- valued  random  noise  vector  Srit  of  length  2N  is  generated.  Then 
the  complex  vector  Sht  of  length  2N  is  generated  according  to  equation  (47). 
The  cumulative  sum  of  6ht  is  then  computed,  resulting  in  another  vector  of 
length  2N,  that  represents  a  complex  random  walk  process.  The  indices  of 
this  vector  are  reversed,  so  that  random-walk  begins  with  the  last  (2A^th) 
element  and  moves  to  earlier  times.  This  vector  is  called  w.  Then,  the 
output  field  Ct  is  computed  for  each  of  the  N  time  increments  t  by  summing 
the  M  +  I  partial  waves  whose  complex  phases  are  given  by  the  values  of 
w.  Thus,  the  first  sample  of  Ct  is  computed  by  starting  at  the  middle  of 
the  vector  w  (i.e.,  and  using  the  previous  M  values.  The  complex 

phase  of  the  first  partial  wave  at  time  increment  t  is  given  by  Wt+N-i  — 
Wt+(N-2i  and  the  last  partial  wave  has  complex  phase  Wt+N-i—Wt+N-i-{M+i)- 
Then,  the  magnitude-squared  and  phase  of  Ct  are  computed.  The  magnitude- 
squared  of  Ct  is  proportional  to  the  intensity  of  the  output  field  at  time  t. 
The  probability  distribution  (  histogram  )  of  the  intensity  is  also  computed. 
Finally,  a  Hanning  (raised-cosine)  windowing  function  is  applied  to  the  time 
series  Ct  and  a  complex  fast-Fourier-transform  (FFT)  is  performed  on  this 
time  series.  The  magnitude  squared  of  the  FFT  is  the  computed  power 
spectrum  of  the  output  field. 


responding  frequency  responses  Fh{f)  and  where  /  is  the  Fourier  frequency. 

Both  impulse  responses  represent  decaying  exponentials,  leading  to  Lorentzian  fre¬ 
quency  responses.  The  frequency  response  of  the  cavity  to  the  multiplicative  noise 
process  is 


Ui-A-,)r;  1  +  i^ 

and  the  response  to  the  input  field  fluctuations  is 


(2.40) 


1 

1 

UJc 


(2.41) 


where  Uc  =  T/t  is  the  rolloff  frequency  of  the  Lorentzian  response,  and  the  Fourier 
frequency  /  is  the  offset  in  Hz  from  the  optical  frequency  Um  ■ 

The  power  spectrum  of  the  output  field  fluctuations  is  comprised  of  two  terms 
due  to  the  multiplicative  and  input  fluctuations.  The  carrier  at  the  cavity  resonance 
frequency  Um  is  shifted  to  DC  due  to  the  discrete-time  analysis,  and  has  amplitude 


im’ 

This  is  equivalent  to  the  output  power  at  the  peak  of  the  resonant  amplifier  trans¬ 
mission  as  derived  previously  in  equation  (2.5). 

The  noise  processes  6ht  and  SRt  give  rise  to  complex- valued  terms  in  the  out¬ 
put  field.  It  is  assumed  that  the  noise  processes  6ht  and  8Rt  are  uncorrelated,  so 
the  power  spectra  due  to  each  process  may  be  added.  The  term  due  to  multiplica¬ 
tive  noise  is  the  product  of  the  multiplicative  noise  complex  power  spectral  density 
Shif)  and  the  magnitude-squared  of  the  frequency  response  Fk{f).  Similarly,  the 
contribution  due  to  the  input  noise  is  the  product  of  the  input  field  noise  complex 


B-42 


power  spectral  density  SR{f)  and  the  magnitude-squared  of  the  frequency  response 
Fh(/)-  Therefore,  the  total  output  field  noise  power  spectrum  due  to  both  input 
field  fluctuations  and  multiplicative  gain  and  phase  fluctuations  is: 


Seif)  =  lW)r  Snif)  +  \FH{f)\^  Snif)  (2.43) 

Substituting  equations  (2.40)  and  (2.41)  for  the  frequency  responses  Fh{f)  and 
Fh(/),  and  using  the  approximation  when  Ko  ~  1  that  T  =  \ln{Ko)\  ~  1  —  Aq, 
yields  the  total  output  field  power  spectrum: 


ScU)  = 


'r^HoR\  1 


Snif)  + 


27'i/oA  \  1 

(l-A^.)rj  1  +  z^ 


rnrj  rr^  i  +  ^ 


Skif) 


Skif)  . 


(2.44) 


This  is  the  general  expression  for  the  output  field  fluctuation  power  spectrum  of  a 
resonant  amplifier  including  input  field  noise  and  multiplicative  noise. 

Consider  now  the  first  term  of  this  expression  due  to  the  input  noise  6R.  For 


low  gain  and/or  low-level  multiplicative  noise,  the  signal-to- noise  ratio  for  Fourier 
frequencies  less  than  Uc  will  be  the  same  as  the  input  signal-to-noise  ratio.  Physically, 
this  implies  that  the  output  field  phase  and  amplitude  exactly  follow  the  input  field 
phase  and  amplitude  fluctuations  for  fluctuation  rates  less  than  the  cavity  bandwidth 
cUc,  which  is  intuitively  expected.  For  fluctuation  rates  faster  than  Uc,  the  input  field 
fluctuations  are  filtered  by  the  cavity,  and  so  are  rolled  off  as  l/p.  As  the  net  gain 
Ko  is  increased,  the  cavity  corner  frequency  Uc  decreases,  and  the  integrated  output 
noise  power  due  to  the  input  fluctuations  decreases,  implying  a  decreased  spectral 
width  of  the  output  field.  The  resonant  amplifier  therefore  produces  a  narrowband 


output  signal  of  increasing  spectral  purity  from  a  broadband  noise  input  as  the  net 
round-trip  gain  Ko  is  increased.  This  is  the  expected  behavior  of  a  resonant  amplifier 
as  obtained  from  the  standard  analysis  leading  to  equation  (2.5),  and  is  the  basis  of 
the  linear  model  of  the  laser  as  a  noise-driven  resonant  amplifier,  to  which  we  will 
return  later. 

Now  consider  the  second  term  of  equation  (2.44)  .  This  term  due  to  multiplicative 
noise  is  new,  and  not  obtained  from  the  standard  analysis  of  the  resonant  amplifier. 
The  multiplicative  noise  term  is  multiplied  by  the  factor  1/(1  —  KoY  compared  to 
the  first  term  due  to  the  input  field  fluctuations.  It  is  seen  that  as  the  net  gain 
Ko  is  increased,  given  fixed  levels  of  input  fluctuations  SR{f)  and  multiplicative 
fluctuations  Shif),  eventually  the  output  field  power  spectrum  Seif)  become 
dominated  by  the  multiplicative  noise  term.  After  this  happens,  the  signal-to-noise 
ratio  for  Fourier  frequencies  less  than  uJc  will  decrease  with  further  increases  in  the 
net  gain. 

Again,  making  use  of  the  fact  that  for  /C  ~  1,  F  =  \in{Ko)\  ~  (1  —  A'o),  the  roll¬ 
off  frequency  may  be  approximated  as  ~  (1  —  Ko)!’’'-  Making  this  substitution 
in  the  second  term  of  equation  (2.44)  yields 


Scif)\Ko- 


1 


SrU)  + 


Shif)  ■  (2-45) 


.  (1  -  A:,)2 (27r/r)2 

In  the  limit  as  Ko  — >  1,  the  corner  frequency  >  0,  and  the  second  term  due  to 
multiplicative  noise  overwhelms  the  first  term.  The  output  power  spectrum  then 
becomes 
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fluctuations  and  increasing  linewidth  make  the  large-angle  regime  something 
to  be  avoided,  so  that  a  prediction  of  the  onset  of  the  large  angle  regime  is 
sufficient  for  most  applications. 

By  the  Wiener-Khinchin  Theorem,  the  power  spectrum  of  a  stationary 
noise  process  is  the  Fourier  transform  of  its  auto-correlation  function  [9] .  For 
the  case  of  small- angle  multiplicative  phase  noise,  the  amplitude  fluctuations 
of  the  output  field  are  small  compared  to  the  phase  fluctuations,  so  the  field 
auto-correlation  function  is  determined  primarily  by  the  phase  fluctuations. 
In  the  large-angle  regime,  large  amplitude  and  phase  fluctuations  occur,  as 
seen  in  the  previous  chapter.  In  this  regime,  the  field  auto-correlation  func¬ 
tion  and  linewidth  are  determined  primarily  by  the  amplitude  fluctuations. 

First,  the  amount  of  linewidth  enhancement  is  calculated  in  the  small- 
angle  regime  from  the  power  spectral  density  of  phase  fluctuations  obtained 
from  the  small-angle  analysis.  The  small-angle  fluctuations  are  found  to  lead 
to  an  output-power-independent  Lorentzian  power  spectrum  for  the  resonant 
amplifier.  The  power  spectrum  in  the  large-angle  regime  broadens  with  in¬ 
creasing  net  gain  as  the  numerical  results  of  the  last  chapter  illustrated.  The 
onset  of  the  large-angle  regime  is  derived  from  the  net  round-trip  gain  and 
the  multiplicative  noise  variance. 
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3  Multiplicative  noise  in  laser  linewidth 


3.1  Analysis  of  effect  of  small-angle  multiplica¬ 
tive  noise  on  semiconductor  resonant  amplifier 
linewidth 

As  discussed  in  the  previous  chapter,  the  multiplicative  phase  fluctuations 
in  a  semiconductor  laser  gain  medium  are  larger  than  the  gain  fluctuations 
by  the  factor  /?.  So,  in  the  small-angle  regime,  the  ampUtude  fluctuations 
of  the  field  are  small  compared  to  the  phase  fluctuations.  This  is  illustrated 
graphically  in  the  complex  plane  in  Figure  4a  of  the  previous  chapter.  Prom 
the  numerical  simulation  results  at  the  onset  of  the  large-angle  regime,  it  is 
seen  that  the  phase  fluctuations  are  larger  than  the  amplitude  fluctuations 
of  the  field.  This  means  that  the  multipficative  noise  contribution  to  the 
output  power  spectrum  Sc{f)  is  primarily  phase  noise  for  a  semiconductor 
laser  amplifier  in  the  small-angle  regime. 

The  output  field  power  spectrum  is  the  Fourier  transform  of  the  field 
auto-correlation  function.  Following  the  development  of  Peterfnann  [10],  the 
output  field  power  spectrum  may  be  derived  in  terms  of  the  frequency  fluctu¬ 
ations  of  the  field  when  the  amplitude  fluctuations  are  neghgible.  In  this  case, 
the  field  autocorrelation  function  may  be  obtained  in  terms  of  the  frequency 
fluctuations  of  the  output  field  [10]  (See  Appendix  for  detailed  derivation): 

{C{t)C*{t  -  T))  =  (48) 

where  P  is  the  average  output  intensity,  {^)  is  the  mean  frequency  offset 
from  the  mode  frequency  Um,  and  the  mean  square  phase  fluctuation  (A0^) 
is  related  to  the  power  spectrmn  of  the  frequency  fluctuations  Sv{f)  as 
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(49) 

The  power  density  spectrum  of  the  field  is  then  obtained  as 

/OO  .  ^ 

{Cit)C*{t  -  T))e-^^^dT.  (50) 

-OO 

These  expressions  relating  the  power  density  spectrum  of  the  field  to  the 
frequency  fluctuation  power  spectrum  are  valid  as  long  as  the  the  phase 
fluctuations  exhibit  a  Gaussian  probability  distribution.  The  derivation  of 
these  results  is  detailed  in  the  Appendix. 

When  the  multiplicative  noise  contribution  to  the  output  power  spec¬ 
trum  Sc{f)  (derived  in  the  previous  chapter)  is  predominantly  phase  noise, 
the  power  spectrum  of  instantaneous  frequency  fluctuations  is  related  to  the 
power  spectrum  of  phase  fluctuations  as  [18] 

ScAf)  =  fSc(f).  (51) 

Substituting  the  previous  result  from  the  small-angle  analysis  for  the  power 
spectrum  of  output  phase  fluctuations,  the  power  spectrum  of  the  output 
frequency  fluctuations  for  net  gain  Ko  approaching  unity  may  be  written 

ScAf)  =  f  (2)(77)5  (“) 

The  factors  of  cancel,  leaving 

ScAf)  =  SM)  (53) 
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The  frequency  fluctuation  power  spectrum  thus  has  the  same  form  as  the 
power  spectrum  of  multiplicative  fluctuations  Sh{f)- 

Consider  the  special  case  when  the  power  spectrum  of  multiplicative  fluc¬ 
tuations  Sfi{f)  Is  white.  Equation  (53)  then  implies  that  the  power  spectrum 
of  output  frequency  fluctuations,  Sc,u{f)i  is  also  white.  In  this  case,  the 
autocorrelation  function  of  the  field  amplitude  is  (See  Appendix) 


{C{t)C*{t  -  T))  =  (54) 

where  tc  =  2/ScA^)-  The  Fourier  transform  of  this  expression  yields  a 
Lorentzian-shaped  spectrum: 


WcA)  = 


2tc{P) 

1  +  ([o’  —  (^)]  tcY 


(55) 


which  is  centered  around  the  optical  frequency  +  A),  and  has  spectral 


width 


,  ,  1  _  5c..(0) 

6u\muU  -  —  - 


Sh{0). 


(27r)^r2 

In  the  small-angle  regime,  the  total  linewidth  has  the  functional  form: 


(56) 


where  the  first  term  on  the  right  is  the  linewidth  obtained  from  the  standard 
analysis  without  multiplicative  noise  [4]. 

In  the  standard  analysis  of  the  resonant  amplifier  without  multiplicative 
noise  [4]  discussed  in  the  previous  chapter,  the  linewidth  is  given  by: 
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For  Ko 


1  —  P1P2GI  1 

P1P2GI  ITT 


1,  this  may  be  written 


(58) 


1-Ko 

TVT 


Some  other  useful  relationships  are: 


(59) 


2r  r  Uc  (1  -  Ko)c 

27rr  TTT  TT  27rnZ 


(60) 


The  total  linewidth  of  the  resonant  amplifer  can  finally  be  expressed  in 
terms  of  the  net  gain  Ko  as  the  sum  of  the  standard  and  multiplicative 
contributions 


oi'\tot  = - 1- 


Si(0) 


(61) 


TTT  (27r)^r^ 

when  ii'o  ~  1-  The  first  term  from  the  standard  analysis  decreases  with 
increasing  net  gain  Ko,  whereas  the  second  term  due  to  multiplicative  noise 
is  approaches  a  constant  as  — >  1. 

A  noise  source  that  has  a  white  power  spectrum  up  to  infinte  frequency 
has  a  delta-function  auto-correlation  function,  which  is  a  useful  mathemati¬ 
cal  model,  but  is  not  physically  realistic.  A  more  realistic  model  for  a  noise 
process  has  an  exponential  time-corellation  proportional  to  with  cor¬ 
relation  constant  a.  The  Fourier  transform  of  an  exponentially-correlated 
noise  source  yields  a  Lorentzian  power  spectral  density,  which  is  white  up  to 
a  Fourier  frequency  of  approximately  f  =  a,  and  then  rolls  off  as  l/P  for 
higher  frequencies.  Then,  for  practical  purposes,  a  white  noise  source  may 
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be  modeled  as  having  a  a  larger  than  all  Fourier  frequencies  of  interest.  In 
practice,  the  noise  may  be  considered  to  be  white  if  a  is  at  least  as  large 
as  the  cavity  mode  spacing  fc,  equation  (56)  may  be  used  to  estimate  the 
linewidth. 

When  a  is  small,  then  the  approximation  of  equation  (56)  is  no  longer 
valid,  and  the  lineshape  of  the  field  power  spectrum  will  not  necessarily  be 
Lorentzian.  The  lineshape  must  then  be  computed  from  the  Fourier  trans¬ 
form  of  the  auto-correlation  function  of  the  electric  field.  The  electric  field 
auto-correlation  function  must  be  first  be  derived  in  terms  of  the  frequency 
fluctuation  power  spectrum,  as  is  done  in  reference  [10]. 

3.2  Onset  of  the  large-angle  regime 

Prom  the  numerical  analysis  of  the  previous  chapter,  it  is  seen  that  when  the 
phase  of  the  later  partial  waves  exceeds  the  small-angle  regime  of  0.1  radians, 
the  output  amplitude  and  phase  fluctuations  increase  dramatically,  (Figures 
2  to  16,  previous  chapter).  In  this  regime,  the  linewidth  of  the  resonant 
amplifier  can  no  longer  detemined  from  the  phase  fluctuations  of  the  output 
field,  as  in  the  small-angle  regime.  As  illustrated  in  Figures  8a-e,  the  power 
spectrum  of  the  electric  field  becomes  progressively  noisier  as  the  net  gain  is 
increased,  corresponding  to  increased  linewidth. 

The  large  amplitude  and  phase  fluctuations  make  this  regime  of  operation 
undesirable,  so  avoiding  the  onset  of  this  regime  is  of  great  practical  interest. 
Prom  the  preceeding  numerical  analysis,  the  onset  of  the  large  amplitude 
fluctuations  just  begins  to  occiu:  for  net  gain  of  Ko  =  0.9  and  when  the 
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multiplicative  noise  standard  deviation  is  c;,  =  0.1  radians,  as  illustrated  in 
Figures  4a  and  5a.  For  this  value  of  net  gain,  the  number  of  partial  waves 
to  achieve  1  percent  convergence  of  the  output  partial  wave  summation  is 
computed  to  be  approximately  45.  So,  the  variance  of  the  phase  for  the  last 
(i.e.,  Mth)  partial  wave  at  the  onset  of  the  large  angle  regime  is 

\/2Mah  «  (0.1)  «  Iradian.  (62) 

So,  to  avoid  rebroadening  of  the  linewidth,  the  net  gain  should  be  kept  below 
the  value  that  causes  the  standard  deviation  of  the  Mth  partial  wave  to 
exceed  the  0.1  radians,  i.e.,  within  the  limits  of  the  small-angle  model,  as 
discussed  in  the  last  chapter. 

3.3  Application  of  multiplicative  noise  model  to  laser 
linewidth 

In  this  section,  the  results  of  the  small  angle  analysis  for  the  effect  of  mul- 
tipUcative  noise  in  a  resonant  amplifier  will  be  applied  to  a  linear  model  of 
a  laser  to  derive  an  analytic  expression  for  the  linewidth.  The  laser  may 
be  considered  to  be  a  resonant  amplifier  driven  by  an  equivalent  input  noise 
source  corresponding  to  the  spontaneous  emission  in  the  active  medium  [19], 
&  [20].  In  this  model,  the  presence  of  the  constant  amplitude  spontaneous 

emission  noise  input  field  causes  the  net  round-trip  gain  to  saturate  at  a 
value  slightly  less  than  unity,  resulting  in  a  finite  linewidth  for  the  resonant 
amplifier.  The  ratio  of  the  output  intensity  to  the  input  noise  intensity  due 
to  spontaneous  emission  determines  the  net  round-trip  gain.  The  wide-band 


spontaneous  emission  noise  is  amplified  within  the  narrow  bandpass  of  the 
resonant  amplifier,  yielding  the  laser  output.  Thus,  for  high  gain,  the  laser 
output  field  may  be  viewed  as  a  wave  of  constant  amplitude  and  frequency,  on 
which  small  statistical  amplitude  and  phase  fluctuations  due  to  the  sponta¬ 
neously  emitted  photons  are  superimposed.  The  modified  Schawlow- Townes 
model  [21]  takes  into  account  the  damping  effect  of  the  saturated  gain  above 
threshold  on  the  amplitude  fluctuations,  which  results  in  a  factor  of  two  de¬ 
crease  in  the  linewidth  above  threshold.  Although  the  linear  laser  model  does 
not  explicitly  include  damping  of  the  amplitude  fluctuations,  its  simplicity 
makes  it  easy  to  see  the  basic  physics  responsible  for  the  linewidth  and  the 
results  for  the  multiplicative  noise  are  easily  obtained. 

In  both  the  original  and  the  modified  Schawlow-Townes  models,  the 
linewidth  at  a  given  output  power  level  is  seen  to  arise  from  the  additive 
noise  source  of  the  spontaneously  emitted  photons  which  add  random  am¬ 
plitude  and  phase  components  to  the  output  field.  The  power  output  of 
the  laser  increases  as  the  saturation  level  of  the  gain  medium  (i.e.,  pump¬ 
ing)  is  increased,  due  to  increased  stimulated  emission.  However,  since  the 
spontaneous  emission  level  is  proportional  to  the  gain,  which  is  essentially 
unity  above  threshold,  so  it  remains  clamped  at  its  saturated  threshold  value. 
Thus,  the  additive  random  phase  contribution  of  the  spontaneous  photons 
becomes  less  significant  as  the  power  increases,  resulting  in  an  increasing 
signal-to-noise  ratio  and  decreasing  linewidth.  Thus,  by  the  the  Schawlow- 
Townes  model,  it  is  predicted  that  the  linewidth  of  a  laser  can  be  made 
arbitrarily  small  by  operating  the  laser  at  an  arbitrarily  high  optical  power 
level.  We  now  derive  the  linear  laser  linewidth  model,  including  the  effect  of 


multiplicative  phase  fluctuations. 

Following  the  development  of  Siegman  [4],  the  output  power  of  the  laser 
P  is  related  to  the  spontaneous  emission  input  noise  power  Pgp  as 

l-RGp  =  2T^  (63) 

where  R  is  the  power  reflection  coefficient  of  the  cavity  mirrors,  so  that 
R  =  T  =  'y'^  =  I  —  R  is  the  power  transmission  coefficient,  and  Gp  is  the 
power  gain  per  round  trip.  The  quantity  1  —  RGp  is  deflned  as  the  “saturation 
amplification”  of  the  resonant  amplifier.  The  spontaneous  emission  power  per 
cavity  mode  can  be  written 

where  h  is  Planck’s  constant,  u  is  the  frequency  of  the  light,  c  is  the  speed 
of  light  in  vacuum,  n  is  the  average  value  of  the  refractive  index,  I  is  the 
length  of  the  cavity,  and  N2  and  Ni  are  the  population  of  the  upper  and 
lower  atomic  levels.  For  low-loss  mirrors,  R^  1,  and  above  lasing  threshold 
the  round-trip  gain  RGp  ^  1,  so  {Gp  -  1)  (1  -  R).  Now,  substituting  this 

and  equation  (64)  into  equation  (63),  and  assuming  lossless  mirrors,  so  that 
T  =  1  —  R,  obtains 


From  the  development  of  the  previous  section,  the  linewidth  of  a  resonant 
amplifier  with  gain  is  given  by  equation  (61) 


The  round-trip  power  gain  RGp  is  approximately  1,  as  is  Kg,  so  we  may 
approximate  1  —  RGp  ~  1  — Kg.  Making  these  substitutions  for  the  saturated 
gain  of  the  laser  into  the  last  expression  then  yields 


c(l  -  R) 
2-Knl 


'  hv  N2  1 

N,  -  Ni  (27r)3T2 


for  the  linewidth  as  a  function  of  the  output  power,  where  6i/c  is  defined  as  the 
“cold-cavity”  linewidth,  or  the  linewidth  of  the  cavity  with  the  gain  medium 
removed.  The  first  term  of  this  expression  is  equivalent  to  .the  original  for¬ 
mula  derived  by  Schawlow  and  Townes  for  the  linewidth  of  the  laser,  which 
is  correct  up  to  threshold,  and  displays  the  familiar  1/P  dependence.  Above 
threshold,  this  term  is  reduced  by  a  factor  of  two  due  to  the  damping  of  the 
amplitude  fluctuations  by  the  saturated  gain  [21].  The  second  term  repre¬ 
sents  the  multiplicative  phase  noise  of  the  laser,  and  is  seen  to  be  independent 
of  the  output  power.  This  equation  predicts  that  the  laser  linewidth  should 
initially  decrease  linearly  with  increasing  output  power,  but  will  finally  reach 
a  minimum  value  proportional  to  the  multiplicative  noise  strength.  If  the 
gain  is  increased  yet  further,  the  linewidth  begin  to  “rebroaden”  at  the  onset 
of  the  large-angle  regime.  A  minimum  linewidth  and  linewidth  rebroadening 
is  typically  observed  in  single- frequency  laser  diodes  at  high  power  levels, 
corresponding  to  high  net  round-trip  gain.  This  model  has  the  same  power- 


independent  character  as  the  linewidth  derived  by  Welford  and  Mooradian 
[11]  using  the  phenomenological  expression  of  equation  (2).  But  the  formal¬ 
ism  used  to  derive  it  also  reveals  that  same  mechanism  can  be  responsible 
for  both  a  power-independent  linewidth  and  linewidth  rebroadening,  as  is 
typically  observed  in  single-mode  semiconductor  lasers. 

It  is  well-known  that  semiconductor  lasers  exhibit  a  power-independent 
linewidth  deviation  from  the  modified  Schawlow-Townes  formula  [11],  [22]. 
As  an  application  of  this  model,  the  minimum  linewidth  of  a  typical  single¬ 
mode  semiconductor  laser  will  be  estimated  due  to  the  effects  of  electron  den¬ 
sity  fluctuations  as  derived  in  the  next  chapter.  Consider  a  single-mode  GaAs 
semiconductor  laser  amplifier  of  length  I  =  250/L/m,  which  implies  round-trip 
time  r  «  5  picoseconds.  Assume  multiplicative  phase  fluctuations  due  to 
electron  density  fluctuations,  as  calculated  in  the  next  chapter,  at  a  level  of 
Sh(j,{0)  =  10"^®  rad^/Hz  for  a  carrier  density  at  threshold  of  10^®  cm“®.  Then 
from  equation  (67),  the  minimum  linewidth  for  the  laser  is  calculated  to  be 
120  Hz.  This  is  30  times  less  than  the  lowest  reported  linewidth  of  3.6  kHz 
for  a  single-mode  laser  diode.  In  other  semiconductor  lasers,  the  minimum 
observed  linewidth  is  2  to  4  orders  of  magnitude  larger.  This  estimate  of  the 
minimum  linewidth  implies  noise  mechanisms  other  than  electron-density 
fluctuations  may  be  dominant  in  the  lasers  tested,  but  that  the  lowest  re¬ 
ported  value  to  date  is  not  far  from  the  minimum  due  to  electron  density 
fluctuations. 

Various  mechanisms  have  been  advanced  to  explain  the  origin  of  the 
power-independent  linewidth  and  linewidth  rebroadening,  including  elec¬ 
tron  number-density  fluctuations  [23],  spatial  hole  burning  [6],  and  spatially- 
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dependent  temperature  and  carrier  fluctuations  [24].  These  or  other  mech¬ 
anisms,  such  as  injection-current  fluctuations,  may  be  operative  in  various 
combinations  in  any  particular  laser.  Whatever  the  source,  multiplicative 
phase  and  gain  fluctuations  in  the  active  medium  will  cause  a  minimum 
linewidth  for  the  output  field  emission  spectrum  due  to  the  fact  that  the  ac¬ 
tive  medium  experiences  thermodynamic  fluctuations  that  alter  the  refractive 
index  [12].  Multiplicative  optical  phase  fluctuations  in  the  active  medium 
due  to  temperature  and  density  fluctuations  arise  from  the  same  consid¬ 
erations  responsible  for  Johnson  noise  in  electronic  components.  Therefore, 
multiplicative  noise  represents  a  fundamental  limit  to  the  minimum  linewidth 
achievable  in  a  laser.  When  this  limit  is  reached,  further  increases  in  gain, 
or  decreases  in  cavity  losses,  will  not  afford  decreased  linewidth.  In  fact,  the 
linewidth  is  predicted  to  rebroaden  if  the  net  gain  is  increased  past  the  point 
that  the  large-angle  regime  is  entered.  Also,  it  is  noted  that  the  total  effect 
of  independent  multiplicative  noise  sources  is  diminished  as  more  sources  are 
added  in  parallel.  This  means  that  the  effects  of  multiplicative  noise  will 
be  more  pronounced  in  structures  with  small  transverse  dimensions,  such  as 
semiconductor  lasers  and  optical  fibers,  as  compared  to  bulk  optical  systems. 

3.4  Conclusion 

An  expression  relating  the  effect  of  the  multiplicative  fluctuations  to  the 
linewidth  of  a  resonant  amplifier  was  derived  using  the  results  of  the  previ¬ 
ous  small-angle  analysis.  The  amount  of  linewidth  enhancement  is  calculated 
in  the  small-angle  regime  from  the  power  spectral  density  of  phase  fluctua- 
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tions  obtained  from  the  small-angle  analysis.  The  small-angle  fluctuations 
are  found  to  lead  to  an  output-power-independent  Lorentzian  power  spec¬ 
trum  for  the  resonant  amplifler.  The  power  spectrum  in  the  large-angle 
regime  broadens  with  increasing  net  gain  as  the  numerical  results  of  the  last 
chapter  illustrated.  The  onset  of  the  large-angle  regime  is  derived  from  the 
net  round-trip  gain  and  the  multiplicative  noise  variance.  For  the  case  of 
multiplicative  phase  noise,  it  is  shown  that  this  effect  produces  a  minimum 
linewidth  of  the  resonant  amplifler,  and  qualitatively  explains  the  phenomena 
of  the  power-indenpendent  linewidth  and  linewidth  rebroadening  in  single¬ 
frequency  semiconductor  lasers. 
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Figure  -  1.  Analysis  Basis  Coi 
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Figure  -  2.  Output  Field  in  Complex  Plane  for  Net  Gain,  Ko  =  0.95 
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Figure  -  9.  Output  Intensity  vs.  Round  Trips,  Ko=  0.99 
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Figure  -  11.  Output  Intensity  vs.  Round  Trips,  Kq=  0.999 
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Figure  -  13.  Ouput  Phase  vs.  Round  Trips,  Ko=  0-95 
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Figure  -  14.  Output  Phase  vs.  Round  Trips,  Ko  =  0.99 
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Figure  -  15.  Output  Phase  vs.  Round  Trips,  Ko  =  0.995 


B  -  74 


500  1000  1500  2000 

t  (round- trios) 


Output  Phase 


Figure  -  16.  Output  Phase  vs.  Round  Trips,  Ko  =  0.999 
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Figure  -  17.  Relative  Probability  vs.  Intensity,  Ko  =  0.90 
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Figure  -  18.  Relative  Probability  vs.  Intensity,  Ko  =  0.95 
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Figure  -  21.  Relative  Probability  vs.  Intensity, 
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Figure  -  22.  Power  Spectrum  vs.  Frequency,  Ko  =  0.90 
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Figure  -  23.  Power  Spectrum  vs.  Frequency,  Ko  =  0.95 
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Figure  -  24.  Power  Spectrum  vs.  Frequency,  Ko  =  0.99 
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Figure  -  25.  Power  Spectrum  vs.  Frequency, 
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MISSION 

OF 

ROME  LABORATORY 


Mission.  The  mission  of  Rome  Laboratory  is  to  advance  the  science  and 
technologies  of  command,  control,  communications  and  intelligence  and  to 
transition  them  into  systems  to  meet  customer  needs.  To  achieve  this, 
Rome  Lab: 


a.  Conducts  vigorous  research,  development  and  test  programs  in  all 
applicable  technologies; 

b.  Transitions  technology  to  current  and  future  systems  to  improve 
operational  capability,  readiness,  and  supportability; 

c.  Provides  a  full  range  of  technical  support  to  Air  Force  Materiel 
Command  product  centers  and  other  Air  Force  organizations; 

d.  Promotes  transfer  of  technology  to  the  private  sector; 

e.  Maintains  leading  edge  technological  expertise  in  the  areas  of 
surveillance,  communications,  command  and  control,  intelligence,  reliability 
science,  electro-magnetic  technology,  photonics,  signal  processing,  and 
computational  science. 


The  thrust  areas  of  technical  competence  include:  Surveillance, 
Communications,  Command  and  Control,  Intelligence,  Signal  Processing, 
Computer  Science  and  Technology,  Electromagnetic  Technology, 
Photonics  and  Reliability  Sciences. 


